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ABSTRACT 


The primary purpose of this contract was to determine the effect 
of biaxial loading on the flaw growth rate of 2219-T87 aluminum 
alloy that would be typical of Space Shuttle cryogenic tankage 
design. In support of this general objective, the stress distri- 
bution and stress concentration factors for several integrally 
stiffened panels under various loading conditions were obtained. 
The flaw growth behavior of both stiffened and unstiffened panels 
under biaxial loading conditions was determined. The effect of 
a complex stress state is studied by introducing flaws in fillet 
areas of biaxially loaded stiffened panels. 
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I. INTRODUCTION 


The principal goal of this work was to determine the effect of 
biaxial loading and stress gradients on the fatigue crack growth 
behavior of 2219-T87 aluminum alloy panels that are typical of 
Space Shuttle tank hardware. In order to satisfy this goal, it 
was necessary to perform experimental and theoretical stress anal 
yses in addition to the experimental fatigue flaw growth behavior 
work. The resulting total program, therefore, consisted of the 
following two-part activity: 

1) Stress Analysis (theoretical and analytical) 

2) Fatigue Flaw Growth Behavior (uniaxial and biaxial) 

To satisfy the program objectives, the following specific tasks 
were performed: 

1) Stress Analysis 

- Theoretical stress analysis of integrally stiffened panels 
using two-dimensional and three-dimensional finite element 
techniques ; 

- Experimental stress analysis of stiffened panels using vari- 
ous loading conditions using photoelastic and strain gage 
techniques ; 

- Experimental stress analysis of unstiffened, welded panels 
to determine welding residual stresses. 

2) Fatigue Flaw Growth Behavior 

Uniaxial flaw growth tests conducted to -provide baseline data; 

Unstiffened and stiffened biaxially loaded panels tested to 
provide flaw growth data; 

Unstiffened, welded panels tested to provide flaw growth data 
under biaxial loading; 

- Comparison of the data with other data using currently avail- 
able flaw propagation analyses. 
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Phase B design studies for the Space Shuttle cryogenic tanks in- 
dicate that the tank walls will have integrally milled stiffener 
construction for both the longitudinal stiffeners and the circum- 
ferential ring stiffeners. The stiffener sizing is based on cri- 
teria that the tank walls will not buckle at a given percent of 
the limit design compression load. The tank walls will also be 
subjected to tension loads at particular times during a mission, 
and the so called "smeared thickness" of the tank walls will be 
based on criteria for carrying limit design tensile loads. 

Extensive experience has been obtained on the design and testing 
of large integrally stiffened tanks from the Saturn and other 
booster programs; however, these tanks were basically designed 
and tested for carrying one-time launch loads and pressures. The 
fatigue behavior for these types of tanks resulting from multiple 
launches and aircraft- type gust and maneuver loads has not been 
adequately investigated. Even though integrally stiffened tanks 
may be reliably made for withstanding one-time ultimate loads 
or burst pressures, the ability to carry thousands of fatigue 
cycles at lower loads and pressures cannot now be confidently 
predicted . 

Fatigue crack propagation testing is the most informative method 
for understanding the effects of different fatigue environments 
on structural designs. Extensive testing has been performed on 
small tensile loaded specimens to obtain basic fatigue crack 
propagation data for 2219-T87 aluminum alloy. Ability to directly 
use this data, however, for analyzing integrally stiffened tank 
walls with complex stress distributions has not been investigated. 

As a result, the empirical approach used in this work was deemed 
best to permit a good engineering assessment of the capability 
of existing analytical methods for prediction of flaw growth be- 
havior. The results of this program should aid in understanding 
the stress distribution in integrally stiffened tankage wall de- 
signs and the fatigue crack growth behavior in these designs. 

In addition, the stress analysis and crack propagation results 
should assist in improving the crack propagation analysis methods, 
and possibly uncover any significant fatigue or crack propagation 
problems that may exist in preliminary Shuttle tankage designs. 
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II. 


PROGRAM PLAN 


The nature of the work performed under each program task is dis 
cussed in this chapter. 


A. STRESS ANALYSIS 


Theoretical - The initial theoretical stress analysis activity 
was used to provide assistance in panel and loading technique 
design. To predict the stress distribution throughout the panel 
and the deflections of the loading straps and grip plates, a 
finite element computer program was used. 

Following finalization of the biaxial panel design, a detailed 
two-dimensional stress analysis was performed to provide displace- 
ments and rotations. The results of the two-dimensional analysis 
were used as inputs for the more refined and precise three-di men- 
sional analysis conducted on two areas of the central bay in 
the vicinity of the stiffeners. 

Experimental -Six stiffened panels (varying skin and stiffener 
thickness combinations) were subject to photoelastic stress anal- 
ysis by Photolastic, Inc. This work performed at our facility in- 
cluded analyses at two stress levels and three loading conditions 
(biaxial ratios). 

Following completion of the photoelastic evaluation, coatings were 
removed and sites for strain gage location were selected after a 
careful review of the photoelastic results. The strain gage evalu- 
ation was performed for two of the stiffened panels. 

Two large biaxial panels containing a long weld through the center 
section were subjected to stress analysis using the technique of 
strain gaging and selective metal removal to relieve residual 
stresses . 
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B. 


FATIGUE FLAW GROWTH BEHAVIOR 


Uniaxial Testing - Uniaxially loaded surface flaw specimens were 
evaluated to provide information for cycles to initiate flaw ex- 
tension from an EDM starter notch, and baseline flaw extension 
data suitable for comparison with other work. Twenty-seven parent 
metal and twelve welded specimens were used for this work. 

Biaxial Testing (Unstiffened Panels) - Cyclic flaw testing was per- 
formed for twenty-two panels under balanced biaxial loading. Each 
specimen, containing four defects was tested with various crack 
shapes, depths, and stress levels. 

Biaxial Testing (Stiffened Panels) - Cyclic flaw testing was per- 
formed for six panels under balanced biaxial loading. Each speci- 
men contained ten defects placed in various locations selected 
on the basis of experimental stress analysis data. 

Biaxial Testing (Welded Panels) - Cyclic flaw testing of weld pan- 
els was conducted at a 2:1 stress ratio to simulate the cylindrical 
portion of a propellant tank. Four panels, each containing four 
defects located in the heat affected zone and weld centerline, 
were evaluated. 
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III. 


MATERIAL AND PROCESSING 


Aluminum alloy sheet and plate (2219-T87) used in this program 
were obtained in the fully heat-treated condition. Four thick- 
nesses [1/8, 1/4, 1/2, and 1 in. (3.2, 6.4, 13, and 25 mm)] of 
stock were used; material for each gage was from a single heat. 

Parent metal panels were machined by mechanical and/or chemical 
milling. The unstiffened biaxial panels were fabricated from 
1/2-in. (13-mm) plate by mechanically milling the reduced section 
to 1/4-in. (6.4-mm) thickness, and then chemically milling to 
a final thickness of 1/16 or 1/8-in. (1.6 or 3.2-mm). This two- 
step operation was required because of difficulties encountered 
during mechanical milling. As material was removed, stiffness 
decreased and the remaining material deflected sufficiently to 
cause irregular cutting and loss of dimensional control. 

Stiffened panels (pockets and grip holes) were mechanically milled 
using numerically controlled equipment. Starting stock was 1— in. 
(25-mm) thick. 

Weld panels for biaxial testing were prepared from 1/2-in. (13-mm) 
stock. Material was chemically milled to 1/4 in. (6.4 mm), welded, 
and then mechanically milled to 1/8 in. (3.2 mm). The resulting 
specimen was identical in appearance to parent metal unstiffened 
panels, having a 1/2-in. (13-mm) thick border for gripping and a 
central reduced thickness section measuring 1/8 in. (3.2 mm), but 
contained a weld extending across the complete width of the speci- 
men. 

Panels were prepared for welding in the following manner. First, 
the aluminum was degreased in trichlorethylene vapor, soaked in 
an alkaline solution for 15 minutes, and deoxidized for 10 minutes. 
Then the edges were filed and after the corners were broken slight- 
ly, the region next to the edge was cleaned with a wire brush. 
Welding was performed longitudinal to the grain direction using 
a direct-current, straight polarity power source and automatic 
equipment. A two-pass butt weld was used. Both passes were made 
from the same side, 2319 filler rod was used for the. second pass. 
The weld schedule used for this work, considered typical for pro- 
duction welding at our facility, is given in Table III-l. 
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Table III-l Weld Schedule, 2219-T87 


Welding Parameter 

1st Pass 

2nd Pass 



Voltage, V 

11.5 

12.5 

Current , amp 

175 

150 

Torch speed, in. /min 

11 

11 

Wire feed, in. /min 

0 

40 

Helium gas coverage, cfh 

70 

70 


Uniaxial test coupons were machined from 1/8-in. (3.2-mm) sheet 
stock for parent metal. For welded coupons, 1/4-in. (6.4-mm) 
stock was welded using the same schedule given above; spec im ens 
were then mechanically milled to 1/8 in. (3.2 mm) after welding. 








IV. 


DESIGN AND STRESS ANALYSIS 


The total effort associated with design, theoretical, and experi- 
mental stress analysis of the biaxial. panels is described in this 
chapter. The design of the unstiffened biaxial panel was performed 
with the aid of finite element analytical techniques to establish 
uniformity of loading and to identify possible areas where malfunc- 
tion might be anticipated during cyclic loading. Design and anal- 
ysis of the attendent fixturing required to apply uniform loads 
was also included in this effort. The stiffened panel design 
was established as identical to the unstiffened design except for 
the incorporation of typical rib stiffening members. The preced- 
ing portion of this work is described in Section A of this chapter. 
The detailed two- and three-dimensional stress analyses of the 
stiffened panel is described in Section B. Experimental stress 
analyses using photoelastic stress analysis and strain gage anal- 
ysis are given in Section C. The experimental analysis to deter- 
mine residual stresses in a welded biaxial panel is also described 
in this section. A comparison of the analytical and experimental 
stress analyses is presented in Section D. 


A. DESIGN 


The biaxial specimen size was selected to provide stiffened regions 
simulating those used in 'actual construction. In order to provide 
a panel containing a rib-stiffened pocket of approximately 8 in. 
square (203 mm square) and introduce loads in a uniform manner, an 
overall size of 28. in. square (711 mm square) was selected. 

The specimen was designed for stressing to 45.0 ksi (311 MN/m 2 ) 
which represents a maximum load of 157.5 kips (705 kN) . Actually, 
the load required was assumed to be closer to 200 kips (896 kN) 
when the stiffened loading region is considered. Based on this 
analysis, the biaxial testing machine was designed to ensure 
200 kips (896 kN) load capability. 

Design and Analysis of Unstiffened Panel 

The panel was analyzed with our in-house finite-element computer 
program. The panel and grip plate were modeled in the computer 
program as flat membranes of different materials and different 
thicknesses, and the load straps were modeled as axial members. 

The results of these computer runs are discussed in detail in the 
following paragraphs. 


\ 
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The basic problem is how to reduce the nonuniform stresses or 
strains in the panel’s center region which arise from three con- 
ditions . 

First, the grip plate experiences a slight amount of bending caus- 
ing less load, and therefore less strain, along the sides than in 
the center of the panel. Second, if the outer two inches of the 
panel are thickened, the resulting strain in this region will be 

even further reduced in the direction parallel to the edge of the 

panel, due to the added stiffness. The third contribution to the 
nonuniform stresses across the center of the panel is that the 
loads are applied at points around the panel and the stresses 
along the panel edge are not uniform or biaxial. This condition 

turns out to produce the least significant effect upon the center 

region because the stresses redistribute themselves very quickly 
within a few inches of the panel edge. However, since the local 
stresses in this outer region are higher than the stresses in the 
center of the panel, fatigue loading could cause this area to fail 
prematurely. 

There are three areas in the panel border region that must be 
watched for critical stresses: (1) the corner where the thin 

inner section meets the thick border sections, an area of high 
shear stress; (2) the bearing stress where the greatest load is 
applied by the loading straps; and (3) the concentrated tensile 
stress around the hole at the midspan of any one side of the pan- 
el. These three stresses have been studied in each model to be 
discussed . 

A detailed description of the six models evaluated to select a 
final configuration is presented in the following paragraphs. 

One-fourth of the test fixture and panel was modeled assuming 
symmetry about two axes. This was done to allow a finer element 
breakdown. Since the actual test fixture is basically symmetrical, 
this assumption is valid. 

The first model considered, shown in Figure IV-la, consisted of 
a panel of constant thickness with no stiffened section around 
the edge of the panel. The load on the straps varied from the 
average load by +1.72% at the center of the panel to -3.45% at 
the outer ends of the panel, a total of 5% variation. The stresses 
along the symmetrical center line of the panel in the direction 
perpendicular to the center line varied from the average stress 
by +0.28% at the center to -0.2% at the edge of the panel, a total 
variation of 0.5%. This would be an ideal panel configuration 
if it were not for bearing stresses. 
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t =0i060 


t =0.125/0.500 



Figure TV-1 Computerized Models > Schematic 
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The second panel configuration, Figure IV-lb, had a thickened 
border section around the edge like a picture frame. This was 
designed to provide sufficient bearing strength and tensile 
strength between holes, as discussed earlier in this report. The 
resulting load variance was +9.48% at the outer end of the panel 
and -5.36% at the center of the panel, a total of 15% due to the 
added stiffness of the picture frame. The stress variation was 
+3.0% at the center of the panel to -5.2% at four inches from the 
edge of the panel. This low stress near the edge was caused by 
the added stiffness of the thicker section along the edge of the 
panel. The peak stress in the corner of the thin section was 
53.5 ksi (369 MN/m 2 ). The peak bearing stress was 54.5 ksi (375 
MN/m 2 ) and the concentrated tensile stress was 60.0 ksi (413 MN/m 2 ), 
using a stress concentration factor of 1.67. The stress would be 
more uniform across the center line of symmetry if we could increase 
the loads at the ends of the panel into the stiffened edges. We 
are limited in the extent to which this may be done by peak bearing 
stresses that are already high, and concentrated tensile stresses 
being high at the edge of the attchment pin holes. It is desir- 
able, however, to have these end loads as large as possible within 
those restraints to have more uniform stresses in the center of the 
panel. This idea tends to discredit the premise that the strap 
loads should be as nearly constant as possible to achieve more con- 
stant stresses in the center region. 

The third and fourth panel configurations represented an attempt 
to use straps bonded along the edges to take out the bearing load 
(Fig. IV-lc and Id). These straps were modeled as thickened ele- 
ments along the edges of the panel. The thickened sections were 
not connected at the corners. On Model No. 3 the loading straps 
were all of the same thickness. The computer run predicted a 
stress of 78.3 ksi (540 MN/m 2 ) at the corner of the panel when 
the corresponding stress at the center of the panel was 42.0 ksi 
(289 MN/m 2 ). In an effort to reduce this stress, the thickness 
of the end loading straps in Model No. 4 were increased by three 
times the thickness of the other loading straps. This was done 
in Model No. 4. The peak stress in the corner was reduced to 
51.3 ksi (353 MN/m 2 ). However, the bearing stress was increased 
from 49.1 ksi (338 MN/m 2 ) to 99.0 ksi (682 MN/m 2 ) for above de- 
sired levels . 

Since the bearing stresses were so critical and it was desirable 
to dump a larger percentage of the load into the edges of the 
panel, we tried a thinner central panel section. Along with this 
change, the integral picture frame border was used to reduce high 
corner shear stresses. This was accomplished in Model No. 5 
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(Fig. IV-le) . The central section of the panel was 0.120 in. 

(3.0 mm) thick, and constant thickness loading straps were used. 

The resulting load variance was -41% at the outer edge to -10% 
at the center. The stress variance was +2.9% at the center to 
-3.1% at four inches from the outer edge. The peak corner stress 
was 52.7 ksi (363 MN/m 2 ), the peak bearing stress was 70.2 ksi 
(483 MN/m 2 ), and the peak concentrated tensile stress was 65.1 ksi 
(448 MN/m 2 ). Comparing these stresses to Model No. 2 shows that 
we have lowered the stress in the corner of the panel by 2.0 ksi 
(13.8 MN/m 2 ), but we have raised the bearing stress by 16.0 ksi 
(110 MN/m 2 ) and the peak tensile stress at the edge of the attach- 
ment holes by 5.0 ksi (34 MN/m 2 ). Using a thinner central section 
appears not to be the correct approach. 

The final run, Model No. 6, had a 0.125-in. (3.2-mm) thick central 
section with 0.5-in. (13'-mm) thick x 2-in. (51-mm) border around 
the edge (Fig. IV-lf ) . The corner was stiffened additionally with 
a thickened triangualar element 1 x 1 in. (25 x 25 mm) . The load- 
ing straps were of constant thickness. The resulting loading 
point variance (L/L avg -1 ) 100 was +14.5% at the ends of the panel 
to -8.1% at the middle of the panel. The stress variance (a/cravg -1 ) 
100 was +5.5% at the center of the panel to -3.4% at four inches 
(102 mm) from the edge of the panel. 

The peak stress in the corner of the panel was 45.9 ksi (316 MN/m 2 ) 
and the peak concentrated tensile stress was 60.8 ksi (418 MN/m 2 ), 
using a stress concentration factor of 1.67. These stress values 
are achieved at a central panel stress of approximately 40 ksi 
(275 MN/m 2 ); balanced biaxial, and are shown in detail in Figure 
IV-2. Although this final run had a high concentrated tensile 
stress, it is very localized and yielding will reduce its value 
after the first cycle. 

There is a tradeoff between the three critical stresses. If the 
bearing stress is increased at the corner of the panel to reduce 
the shear stress in the adjacent thinner corner section, the 
tensile stress around the hole increases. This tensile stress 
around the hole increases if the bearing loads are left alone, 
but the corner of the panel is stiffened to reduce the peak shear 
stress. The increased tensile stress is due to a larger portion 
of the bearing loads near the corner of the panel being carried 
by shear through the stiffened corner section into the thickened 
panel edge. The best choice seems to be the compromise obtained 
in Model No. 6. The three critical stresses — bearing stress at 
the loading holes, net-section tensile stress between holes, and 
the peak stress at the panel corner, appear to be balanced so 
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Figure IV-2 Principal Stress and Load Distribution in Unstiffened 0. 125-in. 
(3.2 mm) Thick Biaxial Quarter-Panel from Finite Element 
Analysis 
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that fatigue failure is equally possible at any one of the three 
locations. Detailed discussion of the situation at each of these 
three locations is given in the following paragraphs. 

The final design for the unstiffened biaxial panel is given in 
Figures IV-3 and IV-4 . 

Design and Analysis of Loading Fixtures 

The loading hardware consists of a plate, straps, and pins. 

Loading Plate - The loading plate provides a load path from the 
single pin joint that attaches to the load train to the 14 pin 
joints that attach the loading straps to the panel. The loading 
plate design is shown in Figure IV-5. The A36 steel plate had 
the same hole configuration as the panel and was designed against 
bearing yield and net section yield using the fatigue endurance 
limit . 

Loading Straps - The 4340 steel loading strap design is shown in 
Figure IV-6. As before, bearing stresses and net section stresses 
were considered. 

Loading Pins - Grade 5 steel bolts were used for pins in order 
to meet the strength requirement for carrying a maximum antici- 
pated fatigue load of 13.7 kips (61 kN) . 

Design of Stiffened Panel 

The stiffened panel is identical with the unstiffened panel except 
for the introduction of the rib stiffener members. The stiffener 
dimensions were arbitrarily selected to be typical of hardware. 
Various combinations of membrane and rib thicknesses were used for 
the six panels. Three rib thicknesses [0.060, 0.125 and 0.250 in. 
(1.5, 3.2, and 6.4 mm)], and two membrane thicknesses [0.060 and 
0.125 in. (1.5 and 3.2 mm)] were used. Figures IV-7 and IV-8 give 
the specifications for the stiffened panel configurations. 


B. THEORETICAL STRESS ANALYSIS 

A finite element stress analysis of the stiffened panel was per- 
formed to provide calculated stresses for comparison with the test 
results and to predict the critical locations for flaws. 
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Figure IV~5 Loading Plate for Biaxial Fatigue Machine 
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The stiffened panel was first modeled in two dimensions with plate 
elements. This step was necessary to reduce the magnitude of the 
problem before making a model in three dimensions with tetrahedral 
elements. Modeling initially with tetrahedrons would have made the 
finite element model too large to handle on the computer if a rea- 
sonable mesh size were used; or if the model were made to fit the 
computer, then the mesh size would be too large to show any mean- 
ingful results. More detailed descriptions of the two-dimensional 
and three-dimensional models follow. 

A finite element model of the stiffened biaxially loaded panel is 
shown in Figure IV-9. Symmetry of the panel made it possible to 
model only one-quarter of the complete assembly. The node point 
mesh was developed to give a large number of elements near the 
center of the panel. The 1/8-in. (3.2 mm) (thin rib) and 1/4-in. 
(6.5 mm) (thick rib) stiffeners were modeled as plate elements 
standing on edge. The seven load fingers that connect to each 
side of the quarter panel are modeled as axial elements linked 
to the fixture arms. Again by symmetry, it was possible to model 
half of one load fixture arm. The fixture model was collapsed 
to a stiffness matrix for the seven node points common to the load 
fingers and one node point for the load application point. The 
collapsed stiffness matrix was then renumbered to fit the load 
finger node points along both edges of the panel model. 

Node points on the X-axis line of symmetry were fully restrained 
against deflection in the Y direction and rotation about the X- 
axis. Node points on the Y-axis line of symmetry were fully re- 
strained against deflection in the X direction and rotation about 
the Y-axis. The stiffeners were fully restrained against rotation 
about the Z-axis at the lines of symmetry. The thin rib was fully 
restrained against deflection in the Y direction at the X-axis; 
the thick rib was fully restrained against deflection in the X 
direction at the Y-axis. The node point at the origin was also 
fully restrained against deflection in the Z direction, thus 
making all other node point deflections relative to the node 
point at the origin. 

Formulation of the elements used in modeling the stiffened panel 
was based on references 1 thru 4. 

All plate elements used in the model have five degrees of freedom 
per node point, or 15 degrees of freedom for each triangular plate 
element, and 20 degrees of freedom for each rectangular plate 
element . 
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A balanced biaxial stress of 30 ksi (207 MN/m 2 ) at the mid-surface 
of the plate was expected from the loading arrangement shown in 
Figure IV-10. The stress analysis showed a biaxial stress of 
about 28 ksi (193 MN/m 2 ) on the stiffener side of the plate and 
32 ksi (220 MN/m 2 ) on the opposite side. The stiffeners, there- 
fore, caused a bending stress of 2 ksi (14 MN/m 2 ) at the center 
of the panel. Biaxial stresses for each element within the central 
bay of the panel balance within ± 1%. Principal stresses were com- 
puted for each plate element; their directions are shown for some 
elements in Figure IV-10. Distribution of the principal stresses 
through the stiffener is shown in Figure IV-11. 

The cross-hatched areas in Figure IV-10 illustrate the locations 
of the segments selected for the more detailed three-dimensional 
analysis using tetrahedral elements. One segment designated the 
"thin-rib" model, cuts across the 1/8-in. (3.2 mm) thick stiffener; 
the other segment designated the "thick-rib" model, cuts across 
the 1/4-in. (6.5 mm) thick stiffener. Each segment was chosen to 
minimize the effect on that segment of the boundary conditions for 
the panel model and the other stiffener. The basis for the boundary 
conditions on the three-dimensional model was the displacement 
(deflection and rotation) calculated in the panel model computer 
run for node points common to the two- and three-dimensional 
models. Deflections for node points on the boundary of the three- 
dimensional model, but not common to the two-dimensional model, 
were calculated by straight-line interpolation. 

Formulation of the tetrahedral element is covered in a report by 
A. Holston, Jr., A Three Dimensional Finite Element 3 Martin Marietta 
Corporation R-71-48637-001, 1971. The element has four node points 
with three translational degrees of freedom at each node point, or 
twelve degrees of freedom for each element. It is a constant 
strain element with thermal strains included. It is essentially 
an extension of the In Plane Triangular Plate Element mentioned 
earlier from two to three dimensions. 

Each segment to be analyzed using tetrahedral elements was sub- 
divided into "bricks" (parallelepipeds) and "wedges." The wedges 
were used in modeling the fillets only. Each brick was composed 
of five tetrahedral elements; wedges were composed of three tetra- 
hedrons (Fig. IV-12) . The thin-rib model had 800 node points and 
2592 tetrahedral elements; the thick-rib model had 750 node points 
and 2432 tetrahedral elements. The material properties for alumi- 
num alloy 2219-T87 (given in Section A) were used. The mesh size 
in the Y direction (thin-rib model) was constant. It was also 
constant in the Z direction for the plate portion of the model, 


IV-16 



) 



BMA 
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- Direction of Max. Tensile Principal Stresses, 

- Areas for Three Dimensional Analysis. 

- Loads in kips. 


Figure IV- 10 Principal Stress Directions and Load Distribution in 0.125-in. 

(3.2 mm) Biaxial Stiffened Quarter Panel from Finite Element 
Analysis 


IV-17 
























































































ELEMENT NUMBER 2324 
185, 190, 258, 259 



ELEMENT NUMBER 597 
185, 189, 190, 258 


Figure IV-12 Tetrahedral Modeling Elements 
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but it varied from fine to coarse in the stiffener. The mesh 
size in the X direction (thin-rib model) varied from coarse at 
the plate edges to fine at the juncture of the plate and stiffener 
(Figure IV-13 and Figure IV-14) . Data for elements in the plate and 
stiffener portion of the model were produced automatically with a 
data generation computer program modified to accommodate the model. 
Element data for one fillet on the thin-rib model were generated 
by hand. These elements were also renumbered to fit the other 
fillet on the thin-rib model and both fillets on the thick-rib 
model. The models were plotted to verify that no errors existed. 
Deflections and rotations for the 14 node points on the mid-sur- 
face of the plate and stiffener common to the two-dimensional 
model and the three-dimensional segments were used to calculate 
the deflections of node points on the cut faces of the segments. 
Rotations of a node point on the mid-surface were converted to 
deflections for node points at a distance from and normal to the 
mid-surface. Deflections for intermediate node points on the cut 
faces of the segments were calculated by straight line interpola- 
tion between the node points common to both two- and three-dimen- 
sional models. 

The two problems (thin-rib and thick-rib models) were solved by 
displacing the boundary node points to the calculated deflections 
and computing the stresses thus induced. Output from the com- 
puter program included six values each of stress and strain, plus 
three principal stresses with their associated principal directions 
for each tetrahedral element in each model. In addition, a value 
known as the distortion energy ratio (DER) was calculated for 
each element. The DER is the ratio of distortion energy in the 
stress state to the distortion energy at yield. It is based on 
the von Mises (J 2 ) theory of plasticity. One computes an effective- 
stress level using the normal and shear stress components and 
then divides by a yield stress. This gives the DER and is a con- 
venient way to compare stress levels among elements. 

An examination of the normal stresses and DERs indicated that ele- 
ments away from the boundaries where node points had been dis- 
placed initially did, in fact, behave independently of the plate 
action enforced at the boundaries. Stresses in the middle tetra- 
hedron and an M outside M tetrahedron for selected bricks were com- 
pared and found to differ only slightly. Figure IV-15 illustrates 
middle and outside tetrahedrons. Normal stresses and DERs were 
nearly symmetric with respect to the rib centerline and also with 
respect to a plane at Y = 1/4 in. (6.4 mm) for the thin-rib model 
(X = -1/4 in. (6.4 mm) for the thick-rib model). A cutting plane 
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was passed through the thin-rib model at Y = 7/32 in. (5.7 mm) 
and the thick-rib model at X = -5/32 in. (4.1 mm) (See Figures 
IV-13 and IV-16). Elements selected as typical for the thin-rib 
model had a Y centroid equal to 7/32-in. (5.7 mm). For the thick- 
rib model, typical elements had a X centroid equal to -5/32 in. 

(4.1 mm). These elements in the thin-rib model also had the same 
X centroids (Y centroids in the thick-rib model) in each "stack" 
of bricks when viewed in the Z direction. Figure IV-12 illustrates 
the elements at the corner of the fillet in the thin-rib model. 

DERs for elements selected from the thin- and thick-rib models are 

tabulated in Figures IV-17 and IV-18. Normal stress (o ) variation 

xx 

in the X direction is plotted in Figure IV-19; normal stress (a ) 

yy 

variation in the X direction is plotted in Figure IV-20. All normal 
and shear stresses for the selected elements are tabulated in 
Appendix A, Table A-l. Further examination of the DERs indicated 
that the critical elements were at the toe of the fillets. These 
are shown as shaded areas in Figures IV-17 and IV-18. 


The stress distribution through the plate thickness at selected 
distances from the stiffener centerline is shown in Figures IV-19 
and IV-20. These may be compared with data theory as follows. 

Two types of classical plate analyses could be performed: one in 

which the stiffeners were neglected and the panel analyzed as an 
isotropic plate; another approach would be to "smear" the stiffeners 
and use orthotropic plate theory. The former would produce con- 
stant stress through the thickness; the latter would show a linear 
variation. Thus, deviations from isotropic and orthotropics plate 
analyses are shown by deviations from constant and linear distribu- 
tions, respectively. Figures IV-19 and IV-20 and Table A-l of 
Appendix A show a peak stress of 41.1 ksi, and the isotropic plate 
theory would give 30 ksi and thus underestimate the peak by 37% 
at the toe of the fillet. Also note that the deviation in o' 


is 


greater than that in a 

yy 


xx 

which is to be expected since this 


stiffener is oriented in the Y direction. 


Figures IV-21 and IV-22 are plots of stresses that occurred in 
elements on the front surface of the thin- and thick-rib segment 
models, respectively. Rear surface plots are given in Figures 
IV- 23 and IV-24. 


In summary, a two-dimensional macroanalysis of the stiffened panel 
was used to establish node point deflections and rotations that 
were used as boundary conditions for a three-dimensional micro- 
analysis of small segments of the stiffened panel. The selected 
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Figure IV- 1? (Sheet 1) Distortion Energy Ratios for Thin-Rib Model 






IV-27 


No. DER No. DER No. DER No. DER No, DER No. DER No. DER No. DER No. DER No. DER 


1289 .019 I 1546 


1285 .019 1542 .019 


1279 .048 1536 .048 


1275 .048 1532 .048 


1269 .073 1526 



Thin Rib - Outboard _ 

No. - Number of Element whose y = 7/32 Plane 
DER - Distortion Energy Ratio 


1265 .073 1522 .073 2572 .086 


1259 .105 1516 


1255 .109 1512 .113 


1249 .207 1506 .209 


2534 .086 

2536 .100 

2509 . 154 

2484 . 149 

2472 .190 2496 .214 2556 .335 


2460 .180 2521 .313 2554 .467 


2584 .565> 


// 


1245 .220 1502 .221*1 


1239 .355 1496 .354 


1235 .349 


1229 .463 


1797 .294 1880 .401 , 1957 AW! 2040 .439 2117 .387 2200 .345 2277 .326 

x x V // / / / / 

1791 .415 1874 .415 1951 .437 '2034 yC483 2111 .480' 2194 .438 2271 .414 

/ //////. 

1787 .409 1870 .416 1947 .434 2030 .481 -2107 /.478 2190 .438 .2267 .414 

' " ' //// 

21 1781 .415 1864 .391 1941 .398 2024 .454 2101 .484' 2184 .421 2261 .403 


£ of Rib 


Figure IV-17 (Sheet 2) (aonel) 











No. DER No. DER No. DER No. DER No. DER No. DER 



2259 .072 

2261 .087 

2224 .092 

2201 .091 

/ / / / 2239 .500 2189 .267 2178 .191 

2278 . 522. 2241 .428 2212 .354 2166 .268 

' / / / 


373 427 .359 510 .349 632 .339 895 .292 

382 421 .377 504 .404 626 .405 889 .419 


No. 

DER 

686 

.004 

682 

.005 

676 

.024 

672 

.024 

666 

.046 

662 

.046 

656 

.080 

652 

.087 

646 

.261 

642 

.237 

636 

.349 

632 

.339 

626 

.405 




No. 

DER 

No. 

DER 

No. 

DER 

1469 

.004 





1465 

.004 





1459 

.024 





1455 

.024 





1449 

.046 





1445 

.046 

2414 

.048 



1439 

.083 

2375 

2377 

.073 

.090 



1435 

.091 

2348 

2325 

.097 
.095 | 



1429 

.269 

2314 

2302 

.198 i 
.272 | 

2337 

2360 

.266 

.343 

1425 

.241 

1637 

.250 

1720 

.274 

1419 

.351 

. 1631 

.358 

1714 

.359 

1415 

.340 

1627 

.348 ; 

1710 

.359 

1409 

.406 

1621 

.405 i 

1704 

.379 



Thick Rib - Outboard __ 

No. - Number of Element whose y = 5/32 Plane 
DER = Distortion Energy Ratio 


/ / / 





IV-28 


























j NORMAL ^ STRESS, 

g.x «•* &« 


°XX' 


Figure IV-19 Plot of Normal Stress (a ) 












Figure IV-21 Surface Element Stresses (Front Face) for Thin-Rib Model 


IV-31 





Lower Surface Element, Normal Stress, ksi 



segments included portions of the thin rib and thick rib. A ratio 
(DER) relating the distortion energy level in a given state of stress 
to the distortion energy level at yield stress was introduced as 
a means of interpreting the voluminous data available from the 
three-dimensional analysis. Based on the DER, it was evident that 
normal stresses were nearly symmetric, but not exactly, in both the 
X and Y directions. Tabulation of the DERs for typical sections 
cut through each model indicated little difference in the magnitude 
of the stresses nor location of the critical elements in the two 
models. The critical stress point in each model was shown to be 
the elements at the toe of the fillets adjacent to the plate por- 
tion of the panel. Plots of the normal stress confirmed the need 
for the three-dimensional analysis around the fillets and the 
juncture of the plate and stiffener. The maximum deviation from 
isotropic plate theory, wherein the stiffeners are neglected, was 
37% at the fillet toe for the conf ieuration shown. This deviation 
diminished to 9% at three stiffener-widths away from stiffener 
centerline. 


EXPERIMENTAL STRESS ANALYSIS 


Three distinct experimental stress analyses were conducted in this 
portion of the program. These were (1) photoelastic stress anal- 
ysis of stiffened biaxial panels; (2) strain gage stress analysis 
of stiffened biaxial panels; (3) strain gage residual stress anal- 
ysis of unstiffened, welded biaxial panels. All three experimental 
analyses were conducted on full size panels. The strain gage and 
photoelastic analyses were performed using the biaxial test fixture 
shown in Figure IV-25 and the stiffened panel (Fig. IV-7) . 

Photoelastic Stress Analysis 

All six stiffened panels were analyzed by the photoelastic tech- 
nique. This work was subcontracted to Photoelastic, Inc., of 
Malvern, Pennsylvania, who are recognized experts in this field. 

A summary of their work is included in this subsection. Detailed 
descriptions of their work and findings are given in Appendix B. 

The purpose of this work was to provide determination of (1) stress 
distribution in the center bay of each panel on both the front and 
rear surfaces; (2) level of stress concentration in the corners 
and fillet radii adjacent to the center bay; (3) stress distribu- 
tion in the stiffener ribs. 
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Figure IV-25 Biaxial Testing Machine with Strain Gaged 3 
Stiffened Panel 




The above information was obtained for each panel for three load- 
ing conditions: 1:0 uniaxial loading, 1:1 biaxial loading and 

2:1 biaxial loading. 

Application of the coating was made by Photolastic at their Malvern 
facility. The coated specimens were returned to Denver and stressed 
in the same equipment used for the cyclic flaw growth work. (A 
description of the testing machine is included in Chapter V.) 

Martin Marietta applied the loads and Photolastic performed the 
stress analysis measurements. 

The following observations were made on the basis of the photo- 
elastic study: 

1) A stress concentration of about 15-30% was noted on the back 
smooth face of the panel directly behind the ribs in a di- 
rection perpendicular to the rib length. This stress concen- 
tration was highest directly behind the point of rib inter- 
section; 

2) The central bay section of all panels is essentially an area 
of plane stress indicating good biaxiality and the absence of 
large bending stresses through the thickness; 

3) Although the corners of intersecting ribs showed steep grad- 
ients, the stresses were much less than the nominally applied 
membrane stresses; 

4) Rib thickness did not appear to influence the magnitude of 
the localized stress disturbance, but did affect the extent 
of the disturbance. 

Results for the photoelastic analysis of panels 1 and 6 have been 
extracted from Table I, Appendix A, and are summarized in Figures 
IV- 2 6 and IV-27 . 

Strain Gage Stress Analysis 

Using the results from the photoelastic analysis, appropriate 
locations and types and sizes of strain gages were selected to 
obtain complementary strain gage data for comparison of these 
two experimental analysis methods. The photoelastic coating was 
removed before strain gage installation. 
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Figure IV-26 Selected Stress Data from Photoelastio and Strain Gage 
Analyses (Panel No. 1, 1:0. Stress Ratio) 
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Figure IV- 27 Selected Stress Data from Photoelastic and Strain Gage 
Analyses ( Panel No. 6j 1:0 Stress Ratio) 
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A prime consideration in gage selection was the stress gradient 
found from the photoelastic results. Large rosette gages were 
used in areas of relatively uniform strain, such as the membrane 
regions of the central bay. Small stacked rosettes were used at 
the fillet junction of the two stiffener ribs to determine princi- 
pal directions. Small uniaxial gage pairs were placed in the 
fillet areas of the stiffener rib and behind the ribs. Uniaxial 
gages were orthogonally oriented. 

Approximately 28 strain gages were applied to each panel. Gage 
sizes varied from 1/4 to 1/32-in. (6.4 to 0.8 mm) grid length. 

All were 120-ohm resistance and were compensated for the thermal 
expansion of aluminum material. A cyanoacrylate room-temperature 
curing adhesive was used to bond all strain gages. Customary 
application procedures were used which include careful surface 
preparation by chemically cleaning and etching before gage instal- 
lation. A typical panel gage installation is shown before test 
in Figure IV-28. Note that all gages are located in the center 
bay of the panel. 

The instrumented stiffened panel was placed in the biaxial test 
fixture for application of uniaxial and biaxial loads. Each panel 
was exercised to the maximum required load level to seat the vari- 
ous linkage and pinhole loading points before data acquisition. 

Using a B&F digital acquisition system, strain gage outputs were 
recorded. The strain gage readings were taken as the panel was 
step loaded to predetermined load levels. (The same loads were 
used in the previous photoelastic study conducted on the same 
panels.) The linearity of readings was thus established by hav- 
ing several strain-versus-applied-load points for each unique 
state of loading. 

Loads were applied to achieve a nominal 30 ksi (207 MN/m 2 ) stress 
level in both uniaxial and balanced biaxial tension. Strain gage 
output data was recorded, reduced, and plotted for each gage ver- 
sus load, A best-fit curve was then drawn through these points 
to eliminate zero shift effects. The strains obtained were then 
used to compute principal stresses using rosette and rectangular 
analysis methods. 

Strain gage data for panels 1 and 6 are presented for 1:0 and 
1:1 load ratios in Figures IV- 2 6 and IV-27, and IV-28, IV-29 and 
IV- 30, respectively. 
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Figure IV- 29 Stress Data from Strain Gage Analysis 
(Panel No . 1:1 Stress Ratio) 
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fStress parallel to section. 


Figure IV- 30 Stress Data from Strain Gage Analysis 
(Panel No . 6j 1:1 Stress Ratio) 
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Analysis of the data shows that an area of high stress concentra- 
tion appears to be on the rib side of the panel in the fillet re- 
gion (location D in the four schematic diagrams cited previously). 
The maximum stress concentration was approximately 40%. A second 
region that exhibited some stress concentration was behind the 
rib on the smooth surface of the panel (location G) . Although 
this region did not show concentration in all cases, a concentra- 
tion of 20% was shown for panel No. 1 under 1:0 loading. 

Biaxial stresses in the central region of the panel (location F) 
were quite consistent with a small variation in stress through the 
thickness of the material of Panel No. 6 under 1:1 loading. 

Residual Stress Measurement 


Residual stresses were determined on two welded biaxial panels. 
These panels were processed in the same manner as those used for 
welded flaw growth data. As noted earlier, these panels were pre- 
pared from 1/2-in. plate that was chem-milled to 1/4 in., welded, 
and then mechanically milled to 1/8 in. 

The technique used for this work was to remove metal surrounding 
miniature strain gages and determine the change in strain attendant 
with the relief of residual stresses. Strain gages with a 1/32-in. 
grid size were bonded in both the longitudinal and transverse di- 
rections (with respect to welding direction) at the weld and heat 
affected zone regions. These gages are sufficiently small to char- 
acterize strain distribution in the region of the weld without 
significantly masking or averaging peak reading, which would occur 
if larger gages were used. 

After application of strain gages to the selected areas and strain 
measurement (Fig. IV-31) , the panel was coated with a chemical 
milling maskant. The maskant was then carefully stripped from 
the region surrounding each gage to allow chemical attack. After 
sufficient material is chem-milled from the region surrounding 
each gage, the process is halted, the coating removed, and the 
strain is again measured. Temperature is carefully monitored at 
the time of each reading.* The chemical milling technique allows 
material to be removed from around the strain gage without intro- 
ducing mechanical stresses or damage to the gages and yet achieves 
the degree of isolation or separation from the surrounding metal 
necessary to obtain accurate residual measurements. 


^Although the strain gages were temperature-compensated for alumi- 
num, the precise nature of the work required additional correction. 


IV-44 


SS-AI 



Experimental data from both panels are summarized in Figure IV-32. 
The data is graphically presented as a function of gage location 
in Figure IV-33. The longitudinal data shows a residual stress 
level of 7.8 ksi (54 M/m 2 ) at the weld centerline. Residual stress 
increases to a peak of over 29 ksi (200 MN/m 2 ) and then starts to 
decrease. These findings are not typical of the results we have 
obtained in other work on aluminum alloys, where the longitudinal 
residual stress is constant (at a peak level of approximately 20 
ksi [138 M/m 2 ]) across the weld zone and into the heat affected 
zone and then starts to decrease. However, in our past work, we 
have measured residual stresses on as-welded material where the 
only material removal has been bead shaving. In this work, re- 
moval of half of the weld thickness would be expected to reduce 
residual stresses. It is interesting to note that although the 
residual stress in the weld centerline was low, the residual 
stress 0.3 in. (7.6 mm) from the centerline (heat affected zone) 
was surprisingly high* A more uniform decrease of residual stress 
would have been anticipated as a result of material removal. 

The longitudinal residual stress along the centerline of the weld 
was uniform in the central portion of the panel, decreasing with 
proximity to the ends of the panel, as expected (Fig. IV-33B) . 

The transverse residual stress was compressive (-11 ksi [-76 M/m 2 ]) 
at the weld centerline, increasing to approximately 20 ksi 
(138 M/m 2 ) tension at a distance of 0.3-in. (7.6 mm) from the 
weld centerline. 

The residual strains measured were used to compute stresses in 
both longitudinal and transverse directions. The effect of the 
biaxial state of stress was considered in these computations. All 
data for the plots in Figure IV-33A were obtained from gages placed 
within the central portion of the weld panel where the longitudinal 
residual stress appeared to be constant. A backside strain reading 
obtained in the center of the panel showed a very small bending 
stress through the thickness which was ignored in computations* 


D. COMPARISON AND DISCUSSION 


A comparison of the theoretical and experimental stress analyses 
showed rather good agreement. Although some of the very sharp 
gradients predicted by the theoretical analysis could not be found 
experimentally, the agreement in some cases (i.e.. Panel No. 1, 

1:0 loading ratio and Panel No. 6, 1:0 loading ratio) was re- 
markably good. 
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Figure IV- 32 Strain Gage Location and Residual Strain Readings 
for Welded 2219-T87 Aluminum Residual Stress 
Measurement 
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The 1:0 loading ratio data are summarized in Figures IV-34 and 
IV 35 for the two experimental techniques. These comparisons 
are restricted to the photoelastic and strain gage methods since 
the theoretical analysis was performed only for biaxial loading. 

Figure IV-34 (Panel No. 1) shows similar behavior in the membrane 
area. . The stress perpendicular to the rib decreases sharply at 
the fillet from both strain gage and photoelastic measurements. 

The stress parallel to the rib shows no variation with respect to 
location according to photoelastic stress measurements; however, 
the strain gage curve shows a decrease in the fillet region. Note 
that the strain gage data are not measured over the same range of 
length as the photoelastic measurements; hence, the basic for com- 
parison is limited. 


Figure IV 35 (Panel No. 6) shows very good agreement for both 
methods with respect to both perpendicular and parallel stresses. 
The agreement for the parallel stresses is within 2 ksi (14 HN/m 2 ) . 

For the 1:1 (balanced biaxial) loading ratio, the comparison is 
restricted to the theoretical finite element analysis and strain 
gage analysis. A comparison with the photoelastic measurements 
cannot be made because this method does not resolve the stress 
values. The photoelastic technique establishes points of high 
stress gradient and principal stress directions; however, stress 
data is given in terms of stress difference (a x - a 2 ) only and 
cannot be used for comparison with the other two methods. 

Data for Panels No. 1 and No. 6 are given in Figures IV-36 and 
respectively. The former panel shows good agreement be- 
tween the theoretical and experimental analyses, particularly 
for the magnitude of the peak stress and the gradient at the be- 
ginning of the fillet. For Panel No. 6, the strain gage analysis 
did not detect the peak at the beginning of the fillet as shown 
by the theoretical analysis. 


In general, the strain gage measurements show lower stress values 
than the theoretical analysis predicted. Several reasons are 
apparent to explain this difference. First, even the smallest 
s ze strain gage, 1/32-in. (0.8 mm) grid, may miss or average sharp 
stress concentrations. Second, the strain gage analysis measures 
surface behavior, whereas the theoretical analysis actually yields 
subsurface data. The centroids of the finite elements closest to 
the surface are 1/64-in. (0.4 mm) from the actual surface. In 
addition, it should be noted that the theoretical two-dimensional 
analysis was performed only for Panel No. 6. 
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Figure IV- 34 Comparison of Surface Stresses (Front Face) by 
Strain Gage and Photoelastic Measurements > 
Panel No. 1 3 1:0 Loading 
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Figure IV- 36 Comparison of Surface Stresses (Front Face) by 
Theoretical and Strain Gage Analyses, Panel 
No. 1, 1:1 Loading 
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Figure IV- 37 Comparison of Surface Stresses (Front Face) 
by Theoretical and Strain Gage Analyses , 
Panel No. 6> 1:1 Loading 
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In summary, for the 1:1 data, we find that the finite element re- 
sults show a sharp stress gradient at the intersection of the fillet 
and membrane. This result was confirmed in Panel No. 1 but not in 
Panel No. 6. The finite element analysis shows the gradient to be 
very narrow; the change in stress in a 1/32-in. (0.8 mm) region is 
12 ksi (83 MN/m 2 ). Hence, it is easy to miss such a sharp gradient 
using experimental techniques. Based on our success in finding 
the gradient in one of our two attempts, we must conclude it does 
occur and the magnitude (approximately 1.4 stress concentration 
factor) is realistic. The gradient may actually be steeper than 
the theoretical analysis predicts, since some of the strain gage 
readings indicate compressive stresses at locations further along 
the fillet. The other region that exhibits stress concentrations 
is the back surface of the panel behind the rib. Both strain gage 
and theoretical predictions show a region of concentration wider 
than in the fillet, but of lower stress magnitude. In this re- 
gion, concentration does not appear to exceed 20%. 

The data from these analyses were used for selection of sites for 
flaw growth measurements made on the six stiffened panels. 

The finite element used is a compatible element; thus, convergence 
of potential energy with decreasing element size is assured. The 
element is a constant stress element, hence its calculated values 
represent average stresses over the element volume. Comparisons 
with other mesh sizes were not made but the agreement with measured 
stress volumes indicates the mesh size was sufficiently fine to 
establish trends for the configuration and loading considered. A 
finer mesh or higher order element would provide a sharper defi- 
nition of stress distribution but computer storage and economic 
restraints must set limitations. 


IV-54 



V. 


FLAW GROWTH BEHAVIOR 


That part of the program dealing with cyclic flaw growth behavior 
for uniaxially and biaxially loaded specimens of both parent metal 
and welded material is presented in this chapter. 


A. EQUIPMENT AND METHODS 


A biaxial testing machine was designed and constructed for eval- 
uation of 28-in. (711 mm) square panels at loads up to 200 kips 
(896 kN). The machine is a simple, welded 20-ft (6.1 m) square 
structure whose principal elements are 24— in. (61 cm) I-beams 
welded together with a 1-in (25 mm) thick gusset plates at each 
comer. Load train details are bolted to this framework. 

The hydraulic system for the machine is closed loop servo control- 
led. Each channel is independently controlled to provide for 
variable loading ratios* 

Alignment of each loading axis was performed by transit measure- 
ments and adjustment of the fixed grip ends. Orthogonality was . 
also confirmed using the transit. Following alignment, an unstif- 
fened panel strain gaged with four rosette and two unidirectional 
gages was installed in the machine and loaded to confirm uniformity 
of stresses in the central portion of the panel and to provide a 
load versus stress calibration. 

The rosette gages were placed at the corners of an imaginary 4-in. 
(102 mm) square in the center of the panel to determine in-plane 
variations. These locations were those where defects would be 
subsequently introduced. The uniaxial gages were located on each 
surface at the panel center to detect bending. 

The panel was incrementally loaded (1:1 stress ratio) in the stress 
range of 20-40 ksi (138-275 M/m 2 ) . The variation in applied load 
averaged less than ±1 percent. The resulting stress field showed 
biaxiality within ±2.3 percent. The principal stress variation 
through the thickness (each side of the panel) was ±2.5 percent. 
This outstanding result is indicative of the precautions taken in 
design, analysis, and fabrication of the machine and specimens. 
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The panel, previously loaded to prove biaxiality, was incre- 
mentally loaded to establish the necessary applied load to give 
a biaxial center panel applied stress of 30 and 40 ksi (207 and 
276 MN/m 2 ) . Stresses were computed from strain data using an 
elastic modulus value of 10.5 x 10 6 psi (721 x 10 3 MN/m 2 ). 

Uniaxial crack growth testing was performed in a 100-kip (448 kN) 
MTS closed loop servo controlled testing machine. Each specimen 
contained two flaws located 3.0 in. (76 mm) apart to ensure 
freedom from interactions. The uniaxial test specimen configura- 
tion is shown in Figure V-l. 

All flaws were introduced using electrodischarge machining (EDM) . 
Because of the criticality of the number of cycles required to 
initiate growth in the biaxial panels, extreme care was taken to 
ensure reproducibility of the starter notches. All work was per- 
formed by a single operator using a single machine. A two-step 
process was used for electrodischarge machining. The initial 
machining was accomplished with a rather blunt tool to provide 
a sufficiently large cavity for fluid circulation and flushing 
action; for machining to final depth, the tool was sharpened to 
a fine tip and reinserted. Figure V-2 shows an example of the 
resulting starter defect. 

It was originally planned to use uniaxial coupons to determine the 
number of cycles to cause crack initiation so that the biaxial 
panels could be tested without prior fatigue sharpening. The 
number of cycles required to cause initiation was subtracted 
from the total number of load cycles applied, and the net dif- 
ference was used for growth rate determinations. As shown by 
the experimental data presented in the Section B, this technique 
was not sufficiently precise; therefore, it was necessary to mark 
the flaw. The technique of varying stress, frequency, and stress 
ratio for marking purposes could not be used for this work. 
Staining was selected as the only available technique. NaOH 
solution (10%) was used as the etchant to provide a subtle stain. 
In order to stain without removing specimens from the testing 
machine, small plastic reservoirs were attached to the flaws; the 
specimens were cycled several times in the presence of the stain- 
ing solution to assure intimate contact with the crack front. 

After removal of the reservoir, the area was cleaned and heated 
until the flaw was dry. Subsequent flaw growth was then readily 
measured. This staining technique has been shown by NASA/MSC 
(Ref 5) and Lockheed (Ref 6) to have no effect on flaw growth 
rate. 
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Figure V-2 
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Macrosection of Electrodischarge Machined Starter 
Defect [one division = 0.001 in. (0.025 rmj] 




B. UNIAXIAL FLAW GROWTH BEHAVIOR 

Uniaxial flaw growth properties were determined for both parent 
metal and welded 2219-T87 material. Detailed data for both con- 
ditions are given in Appendix A, Tables A— 2 and A— 3. 

Parent metal data are presented graphically in Figure V-3. With 
the exception of two data points, scatter was relatively small. 

As noted in the preceding section uniaxial data was originally 
intended to provide information with respect to the number of 
cycles required for flaw initiation. The appropriate data, ex- 
tracted from Appendix A-3, is presented graphically in Figure 
V-4. It is clearly shown that the scatter is too large to permit 
use of this method for predictive purposes. 

Weld data are presented graphically in Figure V-5. In the latter, 
both weld orientations and both defect locations are used to con- 
struct a single curve since the effect of these variables is 
relatively small. 


C. BIAXIAL FLAW GROWTH BEHAVIOR-UNSTIFFENED PANELS 


Biaxial flaw growth properties were determined for both parent 
metal and welded 2219-T87 material. Detailed data are given in 
Appendix A, Tables A-4 and A— 5. 

Parent metal data are presented graphically in Figures V-6 and 
V-7 for the longitudinal and transverse flaw orientations, re- 
spectively. The graphical presentations show least square fits 
for the data (1) plotting all points, and (2) by excluding 
several data points showing large scatter. 

Welded data are presented graphically in Figure V-8. Figures 
V-9 and V-10 show the flaw locations in the welded panels. 
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Figure V-5 Uniaxial Flax Growth Data for Welded Material 
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Figure V-6 Biaxial Flax Growth Data for Unstiffened Parent Metal Panels > 
Longitudinal Flaws ; 1 :1 Stress Ratio 
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Figure V-7 Biaxial Flaw Growth Data for Unstiffened Parent Metal Panels 3 
Transverse Flaws ; 1:1 Stress Ratio 
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Figure V-9 Flaw Locations in Welded Biaxial Fatigue Panels, 

■ Longitudinal Flaws 
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PREFLAW LOCATION DIMENSIONS: Sym about <L's 

Tolerances. 3 ± 0.005 

Figure V-10 Flaw Locations in Welded Biaxial Fatigue Panels > 
Transverse Flaws 
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D. 


CRACK GROWTH BEHAVIOR OF STIFFENED PANELS 


After the six integrally stiffened panels, designed in Chapter IV 
(Fig. IV- 7 ) had been subjected to photoelastic and strain gage 
stress analyses, they were used for flawing and crack growth 
studies. The preceding analyses required that the panels be pre- 
loaded up to a load level corresponding to 30 ksi (207 MN/m 2 ) 
nominal stress in the biaxial test fixture. These nominal loading 
levels were used to develop flaw growth behavior in these same 
panels. All fatigue testing of the stiffened panels was performed 
under balanced biaxial tensile loading at a nominal 30 ksi (207 MN/ 
m 2 ) stress level. This testing was performed at a cyclic stress 
ratio of R = 0 and at a rate of 30 cycles per second. 

Starter flaws were placed in critical stress regions indicated by 
the preceding experimental analyses. These regions include fillet 
areas (both corners and along ribs) and areas in back of the ribs 
on the smooth surface. Flaws were also placed in the ribs them- 
selves, although analyses have indicated this to be a low stress 
region. All flaws were placed within the central bay area of the 
panel. No flaws were placed in the central membrane area because 
an extremely good balanced biaxial field was indicated to be pres- 
ent in this region, and the preceding unstiffened panel tests pro- 
vided sufficient data for such a state of loading. 

Specific areas for flawing are shown in Figures V-ll and V-12. The 
slight modification between the two sets of figures was necessary 
to cover the cases of balanced and unbalanced rib thickness. All 
flaws within a given panel were of the same configuration, i.e., 
a/2C and a/t were constant. Four flaw configurations were used 
in these six panels, 2 a/2C ratios (0.15 and 0.5) and a single 
a/t ratio of 0.5, based on membrane thickness. Ten starter flaws 
were placed in each panel. 

Six of the ten flaws were in regions of relatively constant stress, 
i.e., the stress is not likely to vary greatly across the flaw 
front. These six flaws are represented in sections 4, 5, 7, 8, 9, 
and 10 in Figures V-10 and V-12 and can be identified in subsequent 
tables and figures by these numbers, The growth data for these six 
locations, mentioned above, are given in Table A-6 of Appendix A. 
(See paragraph V-E.4.) 
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SECTION 6-6 


Figure V-ll Flam Location Diagram for Stiffened Panels 1 , 3 and 4 
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The flaws in position 1, 2, 3, and 6 are in fillet regions where 
the stress field is likely to be nonuniform on the basis of our 
stress analysis work* These data are given in Table A- 7 of Ap- 
pendix A. Some data presented in this table are based on measure- 
ments at points on the periphery of the flaw where maximum growth 
occurred* An estimated effective 2C was measured to correspond. 
Other flaws in fillet areas exhibited uniform growth and were 
measured in a normal manner* These areas normally showed lower 
growth rates, indicating subnominal stresses and no stress gra- 
dient * 


Rib 


An example of the type of nonuniform growth exhibited by flaws 
in fillet areas is shown by photograph in Figure V-13* All flaws 
in fillet areas always grew more rapidly in the membrane direction 
at the junction of the membrane and the fillet. This nonuniform 
growth was not obtained from any flaws placed on the smooth back 
face of the panel* 

Membrane 



PANEL 6, POSITION 1 


Figure V-13 Fillet Flaw Showing Nonuniform Growth 
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E. 


DISCUSSION AND DATA COMPARISON 


In this section, we will discuss the following: 

1) Data analysis method; 

2) Parent metal behavior of unstiffened panels; 

3) Weld metal behavior; 

4) Stiffened panel behavior. 

1 . Data Analysis Method 

In the following sections the analysis of the data is based on 
presenting the data as a plot of AK vs da/dN. For this purpose 
the stress intensity range was computed using the equation 

AK = F H Aa^ira/Q [1] 

where F is the front surface correction factor and M is the back 
surface correction factor. Figure V-14 illustrates the crack 
geometry under consideration. 




Figure V-14 Semielliptiaal Surface Flaw in a Plate 
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There have been numerous estimates of the factors F and M by several 
investigators that were reviewed in careful detail by Smith and 
by Shah and Kobayashi in a recent ASME special publication (Ref 7). 
For this work the estimates for F and M presented by Shah and 
Kobayashi were used in the calculations of AK. The detailed com- 
parisons found by Shah and Kobayashi indicate that for most of 
the program, magnification factors compare well with the results 
of Smith (Ref 7) and represent a good estimate for design and 
testing purposes. 

In the computation of AK using Equation [1] for the cases of bi- 
axial loading, it is assumed that the stress component parallel 
to the crack surface has no influence on the stress intensity 
factor. This supposition is correct for cases in which linear 
fracture mechanics is applicable, but is probably not correct in 
a case where the effects of plasticity cannot be ignored. Since 
the stress levels being applied are about 50% of yield and in the 
absence of an analytical method to account for the effect of bi- 
axial stress on AK in the presence of plasticity effects, Equation 
[1] was used in the data analysis. 

Because each test specimen had cracks present with a/2C = 0,15 as 
well as a/2C = 0.50, different amounts of crack growth were ex- 
perienced on the various cracks. An attempt was made to control 
the tests so that the more intense slender cracks did not have 
excessive growth, while trying to maintain an amount of growth 
on the round cracks which was not too small. * In fatiguing enough 
to produce accurately measurable crack growths in the round cracks, 
a comparatively large crack growth was induced in the slender 
cracks. Accordingly, the calculations of AK and da/dN for each 
test were done as follows. 

A value of AK was computed for each crack, using Equation [1] based 
on the dimensions after sharpening which were indicated by stain- 
ing. A second value of AK was computed based on the crack dimen- 
sions at the conclusion of the fatigue test. These two AK values 
were then averaged and used in preparing the AK vs da/dN plots. 

The crack growth rates, da/dN, were approximated by the ratio, 

Aa/AN, where Aa is the measured crack growth and AN is the number 
of cycles for each fatigue test. 

The equations and procedures described here were programmed for 
the computer to reduce the data to the form of AK and Aa/AN, 
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The results of the uniaxial and biaxial fatigue flaw growth tests 
were interpreted in terms of comparison with Forman's equation 
(Ref 8) expressed in the form 

C (AK) n 


da/dN = 


(1-R)K_ - AK 


[ 2 ] 


Wherever possible, a least square fit of Equation [2] to the data 
was performed. For the tests done under this program, the stress 
ratio, R, was zero, so the least square fits were done by rewriting 
Equation [2] in the form 

(Kj - AK) da/dN = C(AK) n 
c 

and taking the natural log of both sides of the equation. This 
produces a problem of least square fitting a straight line in the 
log plane to the data to determine C and n. This was done using 
standard procedures and was programmed for the CDC 6400 computer. 
The following paragraphs present the results of the analyses and 
their interpretations. 

2 . Behavior of Unstiffened Panels (Parent Metal) 

As shown in Section C (Fig. V-6 and V-7), elimination of outlying 
points does affect the shape of the flaw growth curves slightly. 

For the longitudinal direction, elimination of points 1 through 
4 decreases the slope of the curve. Point 5 grew through the plate 
thickness, but its deletion produces no further change in the re- 
sulting 36-point curve. 

In order to least square fit a curve to the longitudinal data, it 
was devided to select a value for corresponding to the speci- 
men orientation. These values were taken from a report by Engstrom 
(Table V-l) . To determine the sensitivity of the curve fit to the 
Kj c value selected, K^ c was varied by plus or minus five percent 

of the nominal value presented in the report. This much variation 
in K produces a shift in the curves on Figure V-6 which is small 

compared to the shift caused by deleting the outlying points. 
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The transverse data was treated in a similar way. The two curves 
shown in Figure V-7 demonstrate the effect of eliminating the out- 
lying points 1 and 2. There is a total 40 data points for the trans- 
verse samples and the curve fitted to all the points is so indicated. 
The lower curve was fitted to the data remaining after deletion of 
points 1 and 2. Point 3 was subsequently deleted because the crack 
grew through the full thickness of the plate and accordingly the 
growth rate value for that crack is uncertain. The absence of this 
point, however, did not modify the 37-point curve any further. 

A further interpretation of these results may be obtained from 
Table V-l which summarizes the results of the least square fits 
to the data in terms of the parameters n and C of the Forman equa- 
tion. It is noted that the greatest change in n occurs as a re- 
sult of deleting the outlying points. The effect of the plus or 
minus five percent shange in K^ is seem to be small on n. It is 

also noted that a small change in n produces a fairly large effect 
on the value of C. In comparing the values for the transverse and 
longitudinal orientations, it is observed that there is not much 
difference in comparable values of n and C. This is born out by 
a comparison of the data of Figures V-6 and V-7 from which it is 
found that the data are indistinguishable from one another. The 
values of n and C found in Table V-l for the two orientations pro- 
duce AK vs da/dN curves that are virtually identical for purposes 
of design in the range of the experiments. Figure V-15 further 
shows the similarity of the two curves. 

It is important to make comparison of these results with uniaxial 
fatigue data in order to assess the effects of biaxial fatigue. 

This may be done in part by reference to Figure V-16 in which AK 
vs da/dN is compared for the two types of loading. It is noted 
that while the data are similar in nature, the uniaxial data does 
fall below the biaxial results; that is, in this range of AK and 
da/dN, a crack subjected to uniaxial fatigue tends to grow faster 
than one in biaxial fatigue by as much as a factor of two or 
three. 

For purposes of overall comparison, Table V-2 presents comparable 
n and C values from different sources to further establish the 
effects of biaxial stress on the fatigue of surface flaws. It is 
noted that the n coefficient is smaller for the biaxial case in- 
dicating the tendency for slower growth. Figure V-17 shows a 
graphical comparison of uniaxial data from three sources. Note 
that Forman's data and the Lockheed data (NAS9-11722) are identical. 
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Table V-l Crack Growth Parameters for Biaxial Flaw Growth 


Defect 

Orientation 

Fracture 
Toughness, K , 
ksi ^in . C 

for indicated number 
of points 

c, 

for indicated number 
of points 

Longitudinal 

44. 2 

(Avg)* 46.2 

48.2 

40 pts 
1.736 
1.780 
1.818 

36 pts 
2.061 
2.105 
2.141 

35 pts 
2.027 
2.049 
2.069 

40 pts 
3.20 
3.06 
2.97 

36 pts 
1.62 
1.56 
1.52 

35 pts 

1.76 

1.77 
1.79 

Transverse 

39.3 

(Avg)* 41.3 

43.3 


38 pts 
1.923 
2.001 
2,060 

37 pts 
2.003 
2.039 
2.070 

40 pts 
2.56 
2.32 
2.17 

38 pts 
1.66 
1.49 
1.38 

37 pts 
1.38 
1.36 
1.35 

40 Points = All Data Points 
For Transverse Specimens: 

37 Pts = All Data - Pts 1 and 

36 Pts = All Data - Pts 1, 2, 

For Longitudinal Specimens: 

36 Pts = All Data - Pts 1, 2, 

35 Pts = All Data - Pts 1, 2, 

Use of n and C values with K_ 

Ic 

in y-in. /cycle 

2 

and 3 

3 , and 4 
3 , 4 , and 5 

and AK units of ksi in. produces da/dN 

of Design Allowable Properties - Fracture of 

\ CR 115388, Prepared by The Boeing Company 
3 k 24, 1972. 

L. Engstrom: Determinatic 
2219-T8? Aluminum Alloy . NASi 
under Contract NAS9-10364, Ta* 


Table V-2 Comparison of Biaxial and Uniaxial Fatigue Data for 2219-T87 
at Room Temperature 



Biaxial 

Uniaxial 




This work 

Boeing 

(Ref. 

Lockheed (Ref ♦ 






NAS9-10364) 


L1722) 


n 

c 

n 

c 

n 

c 

n 

c 

All 

Points 

1.78 

2.32-3.06 

2.34 

1.84 

2.17-2.40 

.71-. 75 

2.5 


Outlying 

Points 









Deleted 

2. 0-2.1 

1.49-1.56 

2.18 

2.36 


* 
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Figure V-15 Comparison of Defect Orientation on 


Flab) Growth Rate for Unstiffened, Biaxial Panels 



STRESS INTENSITY 
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Figure V<-16 Comparison of Stress Ratio on Flaw Growth Rate 




Figure V-17 Comparison of Uniaxial Flaw Growth Rate Data from Various Sources 
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3. Behavior of Welded Panels 


Figure V-18 summarizes AK vs da/dN data for the uniaxial and bi- 
axial fatigue growth of surface cracks in or near a weld in a sheet 
of 2219-T87 aluminum. The data for cracks on the weld centerline 
and in the heat affected zone are all plotted together. In study- 
ing the data, it is difficult to discern any significant distinc- 
tions between the uniaxial and biaxial results for the minimal 
number of data points available. 

In keeping with the work described earlier, an attempt was made to 
least square fit Forman's equation to these data. A difficulty 
arises in doing so because a value of K^ c must be chosen and this 

data is not available because of the extremely tough behavior of 
the weld material. It was decided to select a value of 30 ksi- 
in.h for fitting the fatigue data, remembering that varying K 

does not overwhelmingly change the least square fit curve. Table 
V-3 shows the effect of K T on n and C using the least square 

fitting technique. The curves given in Figures V-5 and V-8 are 
the results of this procedure, and it is noted that the cracks 
subjected to biaxial fatigue tend to grow more slowly until a 
crossover point is reached according to the curves. This is due 
to the higher value of n that was predicted for the biaxial data. 
Since there are not a large number of data points, it is inappro- 
priate to draw any conclusion based on the least square fit curves. 
There appears to be little distinction between the fatigue behav- 
ior of surface cracks in or near welds which depends directly on 
whether or not the loading is biaxial. As usual, the problem is 
complicated by the decreased yield strength of the weld material 
and the material in the heat affected zone, and also the presence 
of residual stresses parallel to the weld bead as well as residual 
stresses that vary in the thickness direction. The theoretical 
stress intensity predictions and stress analyses presently avail- 
able are inadequate to properly cope with these difficulties. 

A comparison of the uniaxial data with Lockheed data from NAS9- 
11722 is given in Figure V-19. 


V-26 



STRESS INTENSITY RANGE, AK, (ksi /In.) 




STRESS INTENSITY RANGE, AK, (ksi /"in. 






Table V-3 Effect of Fracture Toughness Level on Flaw Growth 
Parameters in Welded Material 


Fracture 
Toughness, K , 

Uniaxial Analysis* 

Biaxial Analysis* 

i-C 

ksi \in. 

n 

c 

n 

c 

30 

1.88 

2.41 

3.25 

0.044 

36 

1.98 

2.55 

3.45 

0.038 

41.3 

2.03 

2.73 

3.55 

0.038 

46.2 

2.07 

2.93 

3.61 

0.038 

48.2 

2.08 

3.02 

3.63 

0.039 

^Analysis using least square 

fitting technique 



4. Flaw Growth Behavior of Stiffened Panels 

As described in the preceding Section, the ten flaws placed in 
each stiffened panel were analyzed in two different manners. Flaws 
contained in sections 4, 5, 7, 8, 9, and 10 of Figures V-ll and V-12 
and Table A-6 were assumed to be in regions having a relatively 
constant stress field, which allowed reasonable uniform growth pat- 
terns to develop. Data from these flaws are presented in Figure 
V-20 with the average growth rate best fit curves obtained from 
biaxial and uniaxial tests and for parent material presented pre- 
viously in this section. The points, as plotted, were computed 
using Forman's equation (Ref 8) and experimentally measured stress 
values obtained in the region where each flaw was located. 

Note that in general points are within the two curves plotted for 
comparison from biaxial and uniaxial flaw growth data. Exceptions 
are at the lower end of the curve where growth rates were small and 
highly subject to errors from insufficient flaw growth. Data taken 
from these flaws would be expected to fall between biaxial and un- 
iaxial data results as all points are in regions where stresses are 
either uniaxial or unbalanced biaxial. 

Flaws in positions 1, 2, 3, and 6 (Figs. V-ll and V-12) were lo- 
cated in regions where the stress gradient was expected to be 
changing very rapidly, i.e., where stress concentrations had been 
indicated from experimental analysis. 
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Figure V-20 Flaw Growth Data From Stiffened Panels (Uniform Stress Field) 



The analysis of these data was conducted in the following manner. 
First, a maximum Aa/AN was measured from each flaw. Flaws that 
exhibited uniform growth or growth in a direction normal to the 
surface were also measured. The measured Aa/AN was then used with 
the average growth rate curve in Figure V-20 to obtain an average 
AK. This AK was then used to compute an indicated applied stress. 
Thus, the flaw growth rate is used to compute a maximum stress 
value that can be used to compare with values obtained by stress 
analysis methods previously employed. The resulting stresses are 
then divided by the nominal membrane stress to obtain stress con- 
centration values, K^_, which are presented in Table V-4. 

The table is divided into the four flaw positions that can be re- 
lated from each panel. Flaws at positions 1 and 2 are in similar 
areas, i.e., the corner of the central bay where two fillets inter- 
sect, while flaws at positions 3 and 6 are in the fillet area 
adjacent to ribs. Flaws at positions 1, 3, and 6 are similar in 
that all are perpendicular to the fillet surface at 45° to the 
plane of the panel. 

In general, flaws in positions 3 and 6 show a low stress field 
while flaws in positions 1 and 2 show considerable stress con- 
centration. Stress analysis results indicate that a low stress 
field does occur in the fillet radius which agrees with Table V-4 
results shown at positions 3 and 6. Analysis also shows a sharp 
stress concentration of about 40 percent corresponding very closely 
to flaws in position 2 from the tabulation. 



Table V-4 Stress Concentration Factors Computed from Flaw Growth 

Data and Compared with Experimental and Analytical Values 



Flaw Position 
(Fig. V-ll and V-12) 

K t 

Flaw Growth 

Theoretical 

Experimental 

Panel 2 

1 

1.41 



Panel 6 

1 

0.91 

<0.57* 

0.29 

Panel 5 

1 

1.34 



Panel 1 

2 

0.93 

1.40* 

1.40 

Panel 3 

2 

1.35 



Panel 4 

2 

0.98 



Panel 2 

2 

1.57 



Panel 6 

2 

1.32 

1.40* 

1.29 

Panel 5 

2 

1.45 



Panel 3 

3 

0.85 



Panel 4 

3 

0.77 



Panel 2 

3 

0.75 



Panel 6 

3 

0.88 

0.57 

0.57 

Panel 5 

3 

1.16 



Panel 2 

6 

0.75 



Panel 6 

6 

0.76 

0.67 

0.71 


*The theoretical analysis was not conducted at the corner of the 
panel, represented by flaw positions 1 and 2. Values' shown are 
approximations based on analysis at other fillet locations in the 
panel where both a thin and thick rib was studied. 
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VI. 


SUMMARY AND CONCLUSIONS 


The sample designs used in this program proved to be adequate in 
all respects. No failures of any samples occurred and all samples 
proved to be of adequate size and design to obtain required uni- 
form stress fields of sufficient size in central regions of the 
sample. No delamination of any materials was noted. 

Experimental and theoretical analyses of the stiffened panel con- 
figurations studied yield the following conclusions: 

1) The highest stress concentration occurred on the stiffened 
side of the panel at positions where skin meets the fillet 
radius. The concentration appears to be in the order of 40%. 

2) Another stress concentration, but of smaller magnitude, oc- 
curred on the smooth unstiffened side of the panel in back 
of the ribs but in skin thickness material- The stress con- 
centration factor in this region was about 1.2. 

The photoelastic stress analysis method provided a general picture 
of critical areas to be examined but did not detect the sharp 
stress gradient indicated by analytical and strain gage analysis. 
The uniformity of stress fields and principal stress directions 
were identified by photoelastic analysis. Approximate stress 
levels were obtained by this method in regions of low stress 
gradient. 

Our finite element analysis (three-dimensional) has shown that 
two-dimensional analysis methods would be highly inaccurate for 
predicting critical stresses at geometric anomalies, e.g. fillet 
areas, in an integrally stiffened structure or panel. 

Unfortunately, the analytical finite element analysis was not 
available before the strain gage analysis was conducted. If 
this additional input had been available, it would have aided 
strain gage placement and size selection to an even greater ex- 
tent. The extremely sharp gradient in the fillet area was largely 
unsuspected. Nevertheless, sufficient strain gage results were 
obtained to verify the existence of the strain gradient predicted 
by analytical analysis. 
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A comparison of uniaxial fatigue flaw growth properties (room 
temperature) of 2219-T87 aluminum from several sources with the 
biaxial fatigue flaw growth data obtained in this study reveals 
that the presence of a biaxial state of stress tends to decrease 
flaw growth rate in 2219-T87 aluminum. The decrease is in the 
order of a factor of 2 to 3 in flaw growth rate. 

Flaw growth data from biaxially loaded (2 to 1 stress ratio) 2219- 
T87 weld panels was not sufficient to determine whether this load- 
ing condition causes an increase or decrease from uniaxial loaded 
weldments . 

The flaw growth behavior of stiffened panel configurations, in 
general, confirmed the presence of stress concentrations in those 
critical areas indicated by experimental and analytical analyses, 
namely at the skin-fillet junction. Areas directly behind ribs 
and in the ribs themselves or in the fillet itself were found to 
exhibit slower crack growth and were expected to have lower stress. 

It appears that the current crack growth propagation analysis 
methods satisfactorily predict crack growth ratea in complex 
stress states, such as in fillet areas and rib intersection in 
the stiffened panels studied. 
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APPENDIX A 


TABULATED DATA 




Element 

No. 



k-l-O-' 















Element 

No . 

DER 

Normal Stress 

, ksi 

Shear 

Stress 

ksi 

0 

XX 

a 

yy 

a 

zz 

a 

yz 

a 

zx 

a 

xy 

2385 

0.137 

1.4 

20.3 

1.1 

- 0.8 

1.2 

0.1 

2418 

0.089 

0.7 

15.4 

1.5 

- 1.2 

1.1 

0.1 

2420 

0.098 

0.7 

17.4 

1.9 

- 1.3 

- 0.8 

0.1 

2448 

0.085 

- 0.2 

16.0 

2.6 

- 1.2 

0.5 

0.1 

709 

0.458 

34.5 

35.0 

- 0.4 

- 0.3 

- 0.4 

- 0.2 

715 

0.350 

27.6 

32.0 

- 0.7 

- 0.3 

- 0.4 

- 0.3 

719 

0.356 

27.4 

31.8 

- 1.1 

- 0.3 

- 0.6 

- 0.3 

725 

0.229 

21.6 

28.6 

1.0 

- 0.4 

- 0.7 

- 0.2 

729 

0.220 

22.2 

29.2 

2.2 

- 0.4 

1.5 

- 0.2 

735 

0.111 

0.3 

17.9 

1.6 

- 1.0 

1.4 

0.1 

739 

0.105 

0.9 

18.5 

2.8 

- 1.0 

0.5 

0.1 

745 

0.072 

- 0.1 

15.0 

2.7 

- 1.0 

0.2 

0.1 

749 

0.073 

- 0.2 

14.8 

2.5 

- 1.0 

0.3 

0.1 

755 

0.048 

0.1 

12.1 

1.9 

- 0.9 

0.1 

0 

759 

0.048 

0 

12.0 

1.6 

- 0.9 

0.2 

0 

765 

0.019 

0 

7.3 

0.5 

- 0.8 

0.2 

0 

769 

0.019 

0.1 

7.4 

0.6 

- 0.8 

0.2 

0 

966 

0.460 

34.9 

35.2 

- 0.2 

- 0.3 

0 

0.2 

972 

0.351 

27.8 

32.1 

- 0.6 

- 0.3 

- 0.1 

- 0.3 

976 

0.356 

27.6 

31.9 

- 1.0 

- 0.3 

- 0.3 

- 0.3 

982 

0.228 

21.1 

28.3 

0.7 

- 0.4 

- 0.3 

- 0.2 

986 

0.216 

21.7 

28.9 

2.0 

- 0.4 

0.4 

- 0.2 

992 

0.107 

0.8 

18.2 

1.9 

- 1.0 

0.3 

0.1 

996 

0.104 

1.4 

18.8 

3.0 

- 1.0 

- 0.5 

0.1 

1002 

0.072 

- 0.1 

14.9 

2.7 

- 1.0 

- 0.2 

0.1 

1006 

0.073 

- 0.2 

14.8 

2.5 

- 1.0 

- 0.1 

0.1 

1012 

0.048 

- 0.1 

12.1 

1.9 

- 0.9 

- 0.2 

0 

1016 

0.048 

0 

12.0 

1.6 

- 0.9 

- 0.1 

0 

1022 

0.019 

0 

7.3 

0.5 

- 0.8 

- 0.2 

0 

1026 

0.019 

0.1 

7.4 

0.6 

- 0.8 

- 0.2 

0 

1229 

0.463 

35.0 

35.2 

- 0.2 

- 0.4 

- 0.1 

- 0.3 

1235 

0.349 

27.7 

32.1 

- 0.6 

- 0.4 

- 0.1 

- 0.2 

1239 

0.355 

27.4 

31.8 

- 1.1 

- 0.4 

- 0.2 

- 0.2 

1245 

0.220 

20.3 

28.1 

0.9 

- 0.3 

- 0.4 

- 0.1 

1249 

0.207 

21.0 

28.8 

2.2 

- 0.3 

0 

- 0.1 

1255 

0.109 

0.7 

18.2 

1.6 

- 0.9 

- 0.3 

- 0.1 

1259 

0.105 

1.4 

18.9 

3.0 

- 0.9 

- 0.6 

- 0.1 

1265 

0.073 

- 0.1 

15.0 

2.7 

- 1.0 

- 0.2 

0 

1269 

0.073 

- 0.2 

14.8 

2.4 

- 1.0 

0 

0 

1275 

0.048 

0.1 

12.1 

1.9 

- 0.9 

0.1 

0 

1279 

0.048 

- 0.1 

12.0 

1.6 

- 0.9 

0.1 

0 

1285 

0.019 

0.1 

7.3 

0.5 

- 0.9 

0.2 

0 

1289 

0.019 

0.1 

7.4 

0.6 

- 0.8 

0.2 

0 


Table A-l (oont) 



Element 

No . 

DER 

Normal Stress , ksi 

Shear 

Stress 

, ksi 

a 

XX 

0 

yy 

o 

zz 


a 

zx 

a 

xy 

1486 

0.461 

34.7 

35.1 

- 0.4 

- 0.4 

0.3 

- 0.3 

1492 

0.348 

27.4 

31.9 

- 0.7 

- 0.4 

0.2 

- 0.2 

1496 

0.354 

27.2 

31.7 


- 0.4 

- 0.5 

- 0.2 

1502 

0.221 

20.6 

28.3 


- 0.3 

- 0.7 

- 0.1 

1506 

0.209 

21.3 

29.0 

2.4 

- 0.3 

- 1.1 

- 0.1 

1512 

0.113 

0.2 

17.9 

1.3 

- 0.9 

- 1.4 

- 0.1 

1516 

0.107 

0.9 

18.6 

2.7 

- 0.9 

- 0.5 

- 0.1 

1522 

0.113 

- 0.1 

14.9 

2.7 

- 1.0 

- 0.2 

0 

1526 

0.209 

- 0.2 

14.9 

2.5 

- 1.0 

- 0.3 

0 

1532 

0.221 

- 0.1 

12.1 

1.9 

- 0.9 

- 0.1 

0 

1536 

0.354 

0 

12.0 

1.6 

- 0.9 


0 

1542 

0.348 

0 

7.2 

0.5 

- 0.9 

- 0.2 

0 

1546 

0.461 

0.1 

7.4 

0.6 

- 0.8 

- 0.2 

0 

1704 

0.421 

32.7 

34.8 

0.1 

- 0.5 

0.1 

- 0.5 

1710 

0.403 

27.9 

33.2 

- 2.1 

- 0.5 


- 0.6 

1714 

0.409 

27.8 

33.0 

- 2.5 

- 0.5 

- 1.0 

- 0.6 

1720 

0.286 

17.5 

26.9 

- 3.7 

- 0.8 

- 3.1 

- 1.0 

2460 

0.180 

20.2 

28.1 

3.2 

- 0.4 

- 0.9 

- 0.1 

2472 

0.190 

8.0 

22.1 

- 2.7 

- 0.3 


0 

2484 

0.149 

8.26 

23.3 

0.6 

- 0.4 


0 

2509 

0.154 

1.1 

20.7 

- 0.1 

- 0.5 

IB 

0 

2534 

0.089 

- 0.7 

15.6 

1.4 

- 1.1 

- 1.0 

0 

2536 

0.100 

0.8 

17.6 

2.0 

- 1.1 

0.6 

- 0.1 

2572 

0.086 

- 0.1 

16.1 

2.6 

- 1.1 

- 0.5 

0 

1781 

0.415 

31.7 

34.3 

- 0.4 

- 0.5 

0.2 

- 0.5 

1787 

0.409 

29.2 

33.8 

- 1.4 

- 0.5 

0.1 

- 0.6 

1791 

0.415 

29.0 

33.6 

- 1.9 

- 0.5 

- 0.7 

- 0.6 

1797 

0.294 

23.0 

29.6 

- 1.0 

- 0.8 

- 1.8 

- 1.0 

2496 

0.214 

18.8 

29.6 


- 0.6 

- 3.8 

- 1.1 

2521 

0.313 

24.6 

31.2 


- 0.7 

- 7.0 

- 1.0 

1864 ! 

0.391 

31.2 

33.5 


- 0.2 

0.1 

- 0.6 

1870 

0.416 

31.0 

33.6 ! 

mm 

- 0.2 

- 0.6 

- 0.7 

1874 

0.415 

31.2 

33.8 

- 0.7 

- 0.1 

- 2.1 

- 0.7 

1880 

0.401 

28.6 

29.6 


- 0.5 

- 5.2 

- 1.0 

2554 

0.467 

28.6 

i 

29.6 


- 0.5 

- 9.3 

- 1.0 

| 2556 

0.335 

30.3 

30.1 


- 0.6 

- 8.2 

- 1.0 

1941 

0.398 

i 

32.3 

34.0 


- 0.2 

- 0.8 

- 0.6 

1947 

0.434 

33.6 

34.9 

0.2 

- 0.2 

- 1.7 

- 0.7 

1951 

0.437 

33.8 

35.1 

0.6 

- 0.1 

- 3.2 

- 0.7 

1957 

0.493 

41.0 

35.7 

3.6 

- 0.5 

- 5.8 

- 1.0 

2584 

0.565 

41.0 

35.7 

3.6 

- 0.5 

- 9.9 

- 1.0 
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Table A-l (conol) 


Element 
No . 

DER 

Normal Stress, ksi 

a 

XX 

a 

i n 

a 

zz 

2024 

0.454 

34.0 

33.3 

-1.3 

2030 

0.481 

34.6 

34.5 

-1.3 

2034 

0.483 

34.6 

34.5 

-1.5 

2040 

0.439 

38.7 

34.8 

2.5 

2101 

0.484 

38.1 

35.4 

0.7 

2107 

0.478 

34.4 

34.4 

-1.4 

2111 

0.480 

34.3 

34.3 

-1.6 

2117 

0.387 

31.2 

31.1 

-1.1 

2184 

0.421 

35.0 

33.7 

0.7 

2190 

0.438 

33.5 

33.9 

-0.6 

2194 

0.438 

33.6 

33.9 

-0.7 

2200 

0.345 

31.4 

31.3 

0.9 

2261 

0.403 

32.3 

32.4 

-0.6 

2267 

0.414 

29.6 

31.9 

-2.6 

2271 

0.414 

29.6 

31.9 

-2.6 

2277 

0.326 

28.2 

29.7 

-0.7 


Shear 

Stress 

, ksi 

0 

a 

a 

yz 

zx 

xy 

1 

o 

H* 

-0.5 

-0.4 

0 

-1.8 

-0.6 

0 

-1.9 

-0.6 

1 

O 

H 1 

-1.5 

-0.9 

-0.1 

-0.8 

-0.4 

0 

-1.4 

-0.6 

0 

-1.5 

-0.6 

-0.1 

-0.9 

-0.9 

-0.2 

1.1 

-0.2 

0 

-0.5 

-0.3 

0 

0 

-0.3 

0 

-1.3 

-0.2 

-0.2 

-0.7 

-0.2 

0 

-0.3 

-0.3 

0 

0.2 

-0.3 

0 

0.9 

-0.2 












STAGE I GROWTH 





15- 

>25- 

-30-1 

15- 

•25- 

•30-2 

15- 

•25- 

-40- 

■1 

15- 

•25- 

-40- 

■2 

15- 

•50- 

-30- 

■1 

15- 

•50* 

-30- 

•2 

15- 

•50- 

-40- 

• 1 

15- 

•50. 

-40- 

-2 

15- 

-75- 

-40- 

-1 

15- 

-75- 

-40- 

•2 

15- 

-75- 

-45- 

•1 

15- 

•50- 

-45- 

•2 

50- 

-25- 

-30- 

-1 

50- 

-25- 

-30- 

■2 

50- 

-50- 

-15- 

• 1 

50- 

■50- 

-15- 

•2 

50- 

-50- 

-40- 

•1 

50- 

•50. 

-40- 

•2 

50- 

•75- 

-40- 

-1 

50- 

•75- 

-40- 

•2 

50- 

•75- 

-45- 

•1 


0.0296 

0.2010 

0.0293 

0.2030 

0.0255 

0.2000 

0.0298 

0.1970 

0.0592 

0.4070 

0.0604 

0.3810 

0.0616 

0.4110 

0.0611 

0.0415 

0.0931 

0.5810 

0.0920 

0.6000 

0.0938 

0.5930 

0.0591 

0.4070 

0.0296 

0.0710 

0.0307 

0.0720 

0.0627 

0.1230 

0.0606 

0.1230 

0.0602 

0.1240 

0.0602 

0.1250 

0.0901 

0.1810 

0.0910 

0.1800 

0.0932 

0.1850 


1.13 
1.12 
1.03 12.3 
1.08 13.9 
1.12 13.9 
1.15 13.8 
1.08 19.4 

1.07 19.3 
1.10 26.5 

1.08 26.5 

1.02 35.5 

1.03 21.7 
2.01 6.7 
2.05 

2.40 

2.41 
2.27 12.0 
2.25 12.1 
2.33 15.3 
2.31 14.3 
2.29 16.4 


Kt 

N i 



(cycles) 


STAGE II GROWTH 




CYCLIC GROWTH VALUES 


0.15 1.13 

0.16 1.16 


0.52 | 2.35 
0.50 
0.50 
0.49 2.31 

0.49 2.32 
0.50 2.34 



* First two numbers indicate nominal aspect ratio, a/2C; next two numbers indicate nominal flaw-depth ratio, a/t; next two numbers represent 
nominal applied stress, ksi; two flaws in each sample. 

tMinimum cyclic stress is zero, therefore K <= AK = FMoVira/Q, where F “ 1.0 + 0.12(l-a/2C) 2 , M is the back surface correction factor 
(Forman, R. G.: "Computer Analysis of Two-Dimensional Fatigue Flaw Growth Problem." NASA TMX-58086, Feb 1972); units of K are ksi yin. 


Table A— 2 Flaw Growth Data for Parent MetaZ 0.125-in. 2219-T87 
Aluminum — Uniaxial Loading 









































A- 6 


PU-WT-W-25- 15 
LT 

LB 

PU-WT-H-25-15 


PU-WL-W- 25-30 
LT 


PU-WL-W-50-30 

LT 

LB 

PU-WL-H-25-30 

LT 

LB 

PU-WL-H-50-30 


15 0.86 0.0264 

15 0.86 0.0265 


15 0.52 0.0277 

15 0.52 0.0245 


30 0.92 0.0279 

30 0.92 0.0275 


EDM PREFLAW 


0.1778 0.15 

0.0647 0.41 


0.2063 0.13 

0.0629 0.39 


0.2006 0.14 

0.0643 0.43 


30 0.92 0.0607 

30 0.92 0.0597 


30 0.58 0.0263 

30 0.58 0.0263 


30 0.58 0.0592 

30 0.58 0.0598 


0.4103 0.15 

0.1269 0.47 


0.1909 0.14 

u . 1909 0.14 


0.4073 0.15 

0.1281 0.47 


Q n K+t a. 


1.04 4.6 0.0264 

1.89 3.3 0.0265 


1.11 4.6 0.0277 

1.91 3.1 0.0245 


1.00 9.7 0.0305 

1.95 6.6 0.0299 


1.02 14.8 0.0629 

2.15 9.2 0.0612 


1.10 9.0 0.0269 

1.10 9.0 0.0305 


1.12 13.9 0.0745 

2.24 9.1 0.0612 


STAGE I GROWTH 


0.1809 0.15 1.04 

0.0647 0.41 1.89 

0.2063 0.13 1.11 

0.0629 0.39 1.91 

0.2058 0.15 1.02 

0.P707 0.42 1.93 

0.4128 0.15 1.03 

0.1345 0.46 2.07 

0.2034 0.13 1.09 

0.0738 0.41 2.00 

0.4158 0.18 1.20 

0.1381 0.44 2.13 


bTAGE II GROWTH 


|( cycles)] (in.) 


4.6 420 0.0350 

3.3 420 0.0290 


4.6 400 0.0289 

3.1 400 0.0256 


10.1 2000 0.0523 

6.9 2000 0.0321 


15.0 400 0.1078 

9.5 400 0.0648 


9.2 1000 0.0340 

6.9 1000 0.0325 


15.5 600 0.0958 

9.4 600 0.0653 


V 2C 2 

«2 

Ktt 

X 2 

0.19 

1.15 

5.0 

0.45 

2.08 

3.3 

0.14 

1.12 

4.7 

0.41 

1.98 

3.1 

0.24 

1.27 

11.7 

0.43 

1.98 

7.0 

0.26 

1.31 

20.6 

0.47 

2.15 

9.7 

0.17 

1.17 

9.9 

0.42 

2.02 

7.0 

0.22 : 

1.30 

18.2 

0.45 

2.16 

9.7 


CYCLIC GROWTH VALUES 

a„~a K +K tt A* 

2 2 I j I 2 AN 


|(cycles) | (in.) 


f f ft 


| 0.0041] 9.6 

MM 


| | | | I I I I .1 I I I I I I I I ° Weld 

**Two flaws were placed In each sample: T - Top, B - Boccom: WT - weld transverse to loading direction, WL - weld parallel to loading direction: 
W means flaw placed in centerline at weld, H means flaw placed in heat affected zone: next two numbers in designation represent a/t ratio: last 
two numbers are applied stress in ksi. 

■'•Note yield strength values were: 

17.5 ksi for weld zone, transverse to load 

32.5 ksi for weld zone, parallel to load 

28.7 ksi for heat affected zone, transverse to load 

52.0 ksi for heat affected zone, parallel to load (same as parent metal) 

-•minimum cycle stress is zero, therefore, K - AK = FMoVWQ, where F - 1.0 + 0.12 (1 - a/2C) 2 , M is the back surface correction factor (Forman, 

R. G.: "Computer Analysis of Two-Dimensional Fatigue Flaw Growth Problems." NASA TMX-58086, Feb 1972), units of K are ksiVin. 

*No stain apparent, hence actual rate is greater than measurements indicate. 


Table A-3 Flaw Growth Data for Welded 0.125-in. 2219-T87 
Aluminum — Uniaxial Loading 



STAGE I GROWTH 




STAGE II GROWTH 





CYCLIC GROWTH VALUES 


060-25-40-1 


TLT 

40 

BLL 

40 

TRT 

40 

BRL 

40 

060-25-40-2 


TIL 

40 

BLT 

40 

BRT 

40 

TRL 

40 

060-35-30-1 


TLL 

30 

BLT 

30 

BRT 

30 

TRL 

30 

060-35-30-2 


BLL 

30 

TLT 

30 

BRL 

30 

TRT 

30 

060-35-40-1 


BLT 

40 

TLL 

40 

BRT 

40 

TRL 

40 

060-35-40-2 


BLL 

40 

TLT 

40 

BRL 

40 

TRT 

40 

060-50-30-1 


TLL 

30 

BLT 

30 

TRL 

30 

BRT 

30 


0.44 

2.04 

0.41 

1.94 


9.4 0.0180 0.1004 

9.2 0.0180 0.0990 

6.7 0.0162 0.0400 

6.6 0.0160 0.0395 

9.1 0.0225 0.0994 

9.0 0.0184 0.0962 

6.1 0.0184 0.0398 

6.3 0.0170 0.0406 

0.0200 0.1474 
No Growth 
No Growth | 

No Growth 

7.5 0.0240 0.1450 

7.9 0.0250 0.1575 

5.5 0.0215 0.0508 

5.4 0.0220 0.0480 

10.7 0.0260 0.1495 

10.9 0.0235 0.1510 

7.4 0.0225 0.0500 

7.5 0.0215 0.0515 

10.7 0.0290 0.1430 

9.3 0.0175 0.1205 

7.1 0.0233 0.0380 
7.2 : 0.0216 0.0472 

9.9 0.0426 0.2120 

9.4 0.0393 0.2030 

6.4 0.0319 0.0686 
6.4 0.0315 0.0665 



760 > 0.0390 0.1200 

760 0.0435 0.1357 

760 0.0215 0.0480 

760 0.0203 0.0467 

3550 0.0457 0.1330 

3550 0.0380 0.1144 

3550 0.0234 0.0490 

3550 0.0235 0.0514 

350 0.0504 0.1730 

350 0.0475 0.1620 

350 0.0232 0.0470 

350 0.0244 0.0600 

4600 0.0470 0.1575 

4600 0.0425 0.1605 

4600 0.0245 0.0535 

4600 0.0240 0.0505 

1100 0.0323 0.1628 
1100 0.0358 0.1595 

1100 0.0252 0.0540 

1100 0.0235 0.0540 

1350 0.0380 0.1506 

1350 0.0220 0.1220 

1350 0.0252 0.0520 

1350 0.0230 0.0520 

2100 0.0605* 0.2600 

2100 0.0580 0.2200 

2100 0.0360 0.0760 
2100 0.0352 0.0760 


1.66 13.5 

1.62 11.9 


1.55 
1.44 
2.20 
2.28 5.7 

1.39 15.1 

1.27 13.7 

2.18 7.9 

2.04 


tFirst two letters orient flaw in relation to sample number on panel (T - top, B *» bottom, L = left, R ■ right); last letter shows flaw 
orientation with respect to grain direction (T *> transverse, L = longitudinal). 

t+First number represents material thickness in thousandths of an inch; secon d num ber represents nominal a/t ratio; third number indicates 
stress level in ksi. Minimum cyclic stress is zero, therefore, K = AK = FMo yrra/Q , where F » 1.0 + 0.12(l-a/2C) 2 , M is the back surface 
correction factor (Forman, R. G.; "Computer Analysis of Two-Dimensional Fatigue Flaw Growth Problem." NASA TMX-58086, Feb 1972); units 
of K are ksi^uT. 

♦Indicates thru flaw. 



0.0210 11.1 

0.0255 11.8 

0.0053 7.2 

0.0043 7.1 

0.0232 11.9 

0.0196 10.8 

0.0050 7.1 

0.0065 



0.0174 13.5 
0.0123 12.6 
0.0027 7.8 
0.0020 7.9 

0.0090 13.2 
0.0045 10.4 
0.0019 7.6 
0.0014 7.5 

0.0179 15.3 
0.0187 13.3 
0.0041 6.8 
0.0037 6.8 



Table A-4 Fatigue Flaw Growth Data for Parent Metal 0. 125-in. 
2219-T87 Aluminum-Balanced Biaxial Loading , Un- 
stiffened Panels 



































A- 8 






EDM PREFLAW 




STAGE I GROWTH 




STAGE 11 

GROWTH 



CYCLIC GROWTH 

VALUES 

Paneli 

Number 

0 

a/a 

ys 

a 

0 


m 

% 

19 


2c i 

a l /2C l 


it ft 

1 

N i 

a 2 

2C 2 

a 2 /2C 2 

a 

11 



ma 

IjlM 

(ksi) 


(in.) 

(in.) 




(in.) 

(in.) 




(cycles) 

(in. ) 

(in.) 




[cycles) 

(in.) 

2 


060-50-30-2 























BLL 

30 

0.58 

0.0280 

0.2130 

0.13 

1.09 

9.7 

0.0370 

0.2206 

0.17 

1.17 

11.3 

2400 

0.0515 

0 . 2344 

0.22 

1.31 

14.8 

1650 

0.0145 

13.0 

8.8 

TLT 

30 

0.58 

0.0275 

0.2190 

0.13 

1.08 

9.7 

0.0334 

0.2233 

0.15 

1.13 

10.6 

2400 

0.0478 

0.2325 

0.21 

1.27 

13.4 

1650 

0.0144 

12.0 

8.7 

BRL 

30 

0.58 

0.0278 

0.0625 

0.44 

2.14 

6.3 

0.0295 

0.0658 

0.45 

2.15 

6.5 

2400 

0.0316 

0.0704 

0.45 

2.15 

6.8 

1650 

0.0021 

6.6 

1.3 

TRT 

30 

0.58 

0.0260 

0.0610 

0.43 

2.05 

6.3 

0.0270 

0.0635 

0.43 

2.05 

6.4 

2400 

0.0290 

0.0660 

0.44 

2.11 

6.5 

1650 

0.0020 

6.5 

1.2 

060-50-40-1 























TLL 

40 

0.77 

0.0301 

0.2055 

0.15 

1.07 

13.5 

0.0426 

0.2152 

0.20 

1.19 

16.0 

1050 

0.0565 

0.2351 

0.24 

1.31 

20.2 

750 

0.0139 

18.1 

18.5 

BLT 

40 

0.77 

0.0283 

0.2320 

0.13 

1.04 

13.2 

0.0426 

0.2055 

0.21 

1.22 

15.6 

1050 

0.0582 

0.2492 

0.23 

1.29 

20.8 

750 

0.0156 

18.2 

20.8 

BRT 

40 

0.77 

0.0300 

0.0695 

0.43 

2.02 

9.0 

0.0320 

0.0730 

0.48 

2.05 

9.3 

1050 

0.0374 

0.0822 

0.46 

2.13 

10.0 

750 

0.0054 

9.6 

7.2 

TRL 

40 

0.77 

0.0295 

0.0640 

0.46 

2.15 

8.6 

0.0320 

0.0698 

0.46 

2.14 

9.1 

1050 

0.0350 

0.0766 

0.46 

2.14 

9.6 

750 

0.0030 

9.3 

4.0 

060-50-40-2 

BLL 

40 

0.77 

0.0295 

0.2175 

0.14 

1.05 

13.5 

0.0314 

0.2175 

0.14 

1.06 

13.9 

500 

0.0474 

0.2292 

0.21 

1.22 

17.2 

1000 

0.0160 

15.5 

16.0 

TLT 

40 

0.77 

0.0280 

0.2170 

0.13 

1.03 

13.2 

0.0298 

0.2200 

0.14 

1.05 

13.6 

500 

0.0450 

0.2300 

0.20 

1.19 

16.9 

1000 

0.0152 

15.3 

15.2 

BRL 

40 

0.77 

0.0284 

0.0653 

0.43 

2.04 

8.8 

0.0295 

0.0670 

0.44 

‘2.06 

8.9 

500 

0.0318 

0.0720 

0.44 

2.07 

9.2 

1000 

0.0023 

9.0 

2.3 

TRT 

40 

0.77 

0.0290 

0.0660 

0.44 

2.06 

8.8 

0.0308 

0.0680 

0.45 

2.12 

8.9 

500 

0.0330 

0.0740 

0.45 

2.09 

9.4 

1000 

0.0022 

9.2 

2.2 

125-25-30-1 























TLT 

30 

0.58 

0.0278 

0.2186 

0.13 

1.08 

9.4 

0.0316 

0.2285 

0.14 

1.11 

9.9 

4400 

0.0650 

0.2345 

0.28 

1.48 

12.2 

1640 

0.0334 

11.0 

20.4 

TRL 

30 

0.58 

0.0274 

0.2092 

0.13 

1.09 

9.3 

0.0310 

0.2161 

0.14 

1.12 

9.7 

4400 

0.0615 

0.2220 

0.28 

1.48 

11.8 

1640 

0.0305 

10.8 

18.6 

BLT 

30 

0.58 

0.0284 

0.0620 

0.46 

2.20 

6.3 

0.0323 

0.0715 

0.45 

2.17 

6.7 

4400 

0.0364 

0.6765 

0.05 

0.97 

11.7 

1640 

0.0041 

9.2 

2.5 

BRL 

30 

0.58 

0.0285 

0.0625 

0.46 

2.19 

6.3 

0.0324 

0.0705 

0.46 

2.20 

6.7 

4400 

0.0384 

0.0785 

0.49 

2.34 

7.0 

1640 

0.0060 

6.9 

3.7 

125-25-30-2 

TRL 

30 

0.58 

0.0320 

0.2210 

0.14 

1.12 

9.9 

0.0580 

0.2520 

0.23 

1.34 

12.1 

6250 

0.0940 

0.2900 

.0.32 

1.64 

14.8 

1000 

0.0360 

13.5 

36.0 

TLT 

30 

0.58 

0.0275 

0.2062 

0.13 

1.10 

9.2 

0.0335 

0.2128 

0.16 

1.15 

9.9 

6250 

0.0458 

0.2150 

0.21 

1.29 

10.9 

1000 

0.0123 

10.4 

12.3 

BRL 

30 

0.58 

0.0298 

0.0697 

0.43 

2.06 

6.7 

0.0343 

0.0795 

0.43 

2.08 

7.1 

6250 

0.0394 

0.0850 

0.46 

2.22 

7.3 

1000 

0.0051 

7.2 

5.1 

BLT 

30 

0.58 

0.0278 

0.0674 

0.41 

2.00 

6.5 

0.0314 

0.0750 

0.42 

2.02 

6.9 

6250 

0.0360 

0.0806 

0.47 

2.14 

7.2 

1000 

0.0046 

7.0 

4.6 

125-25-40-1 























(No data 























Sample not 
Stained) 

125-25-40-2 























BLT 

40 

0.77 

0.0271 

0.2048 

0.13 

1.04 

12.6 

0.0330 

0.2092 

0.16 

1.09 

13.5 

1030 

0.0462 

0.2104 

0.22 

1.25 

14.8 

720 

0.0132 

14.1 

18.3 

TLL 

40 

0.77 

0.0311 

0.2140 

0.15 

1.07 

13.3 

0.0340 

0.2185 

0.16 

1.10 

13.7 

1030 

0.0505 

0.2220 

0.23 

1.27 

15.4 

720 

0.0165 

14.5 

22.9 

TRL 

40 

0.77 

0.0315 

0.0660 

0.48 

2.23 

8.7 

0.0327 

0.0680 

0.48 

2.25 

8.9 

1030 

0.0350 

0.0716 

0.49 

2.29 

9.1 

720 

0.0023 

8.9 

3.2 

BRT 

40 

0.77 

0.0300 

0.0640 

0.47' 

2.19 

8.6 

0.0325 

0.0684 

0.48 

2.22 

8.9 

1030 

0.0346 

0.0735 

0.47 

2.20 

9.2 

720 

0.0021 

9.0 

2.9 


tFirst two letters orient flaw in relation to sample number on panel (T = top, B = bottom, L = left, R = right); last letter shows flaw 
orientation with respect to grain direction (T - transverse, L = longitudinal). 

ttFirst number represents material thickness in thousandths of an inch; secon d num ber represents nominal a/t ratio; third number indicates 
stress level in ksi. Minimum cyclic stress is zero, therefore, K = AK = FMc Vna/Q, where F = 1.0 + 0 . 12 (l-a/2C) 2 , M is the back surface 
correction factor (Forman, R. G.: ,: Computer Analysis of Two-Dimensional Fatigue Flaw Growth Problem. n NASA TMX-58086, Feb 1972); units 

of K are ksi V*n^ 

^Indicates thru flaw. 


Table A- 4 (oont) 











EDM PRE 

FLAW 




STAGE I GROWTH 




STAGE II 

GROWTH 


CYCLIC GROWTH 

VALUES 

Panel'** 

Number 

0 

c/a 

ys 

a 

0 

2C „ 

a/2C 
1 0 


v + 

0 



a l /2C l 

H 



a 2 

2C 2 


119 


N 2 orAN 


I88jg| 

4a 

(ksi) 


(in.) 

(in.) 







1 


P| 

(in.) 

(in.) 

1 

1 


(cycles) 


mm 

\cyclef 

125-35-30-2 




















MjBSi 


■■ 

TLT 

30 

0.58 

0.0465 

0.2905 

0.16 

1.15 

11.8 

0.0490 

0.2957 

0.17 

1.17 

12.0 

2700 ' 

0.0770 

0.3072 

0.25 

1.40 


1600 ; 

0.0280 

13.1 


BLL 

30 

0.58 

0.0410 

0.2904 

0.14 

1.11 

11.3 

0.0490 

0.2935 

0.17 

1.17 

12.0 

2700 

0.0810 

0.3075 

0.26 

1.44 

|rjf| 

1600 

0.0320 

13.2 

KSB 

TRT 

30 

0.58 

0.0390 

0.0980 

0.40 

1.93 

7.9 

0.0413 

0.1025 

0.40 

1.95 

8.1 

2700 - 

0.0482 

0.1100 

0.44 

2.11 


1600 

0.0069 

8.2 

wsm 

BRL 

30 

0.58 

0.0382 

0.0940 

0.41 

1.97 

7.7 

0.0402 

0.1000 

0.40 

1.95 

8.0 

2700 

0.0468 

0.1085 

0.43 

2.08 


1600 1 

0.0066 

8.2 

msm 

125-35-40-1 























TLT 

40 

0.77 

0.0395 

0.2750 

0.14 

1.06 

15.1 

0.0420 

0.2802 

0.15 

1.08 

15.5 

1050 

0.0692 

0.2874 

0.24 

1.31 

18.2 

800 

0.0272 

16.9 

34.0 

BLL 

40 

0.77 . 

0.0400 

0.2775 

0.14 

1.06 

15.2 

0.0525 

0.2960 

0.18 

1.14 

16.9 

1050 , 

0.0870 

0.3280 

0.27 

1.39 

21.0 

800 ‘ 

0.0345 

18.9 

43.1 

TRT 

40 

0.77 

0.0382 

0.0972 

0.39 

1.86 

10.7 

0.0420 

0. 1030 

0.41 

1.92 

11.0 

1050 

0.0490 

0.1145 

0.43 

2.01 

11.6 

800 

0.0070 

11.3 

8.8 

BRL 

40 

0.77 

0.0382 


0.38 

1.80 

10.8 

0.0410 

0.1060 

0.39 

1.83 

11.1 

1050 

0.0495 

0.1175 

0.42 

1.98 

11.7 

800 

0.0085 

11.4 

10.6 

125-50-30-1 























BLL 

30 

0.58 

0.0616 

0.4161 

0.15 

1.13 

14.1 

0.0900 

0.4255 

0.21 

1.28 

17.0 

1600 

0.0990 

0.4416 

0.22 

1.32 

18.4 

400 

0.0090 

17.7 

22.5 

TLT 

30 

0.58 

0.0710 

0.4355 

0.16 

1.16 

15.3 

0.1010 

0.4560 

0.22 

1.31 

19.4 

1600 

0.1103 

0.4786 

0.23 

1.34 

21.6 

400 j 

0.0093 

20.5 

23.3 

TRT 

30 

0.58 

0.0627 

0.1358 

0.46 

2.21 

9.3 

0.0653 

0.1442 

0.45 

2.17 

9.7 

1600 

0.0692 

0.1521 

0.46 

2.18 

10.0 

400 

0.0039 

9.8 

9.8 

BRL 

30 

0.58 

0.0620 

0.1335 

0.46 

2.23 

9.3 

0.0647 


0.46 

2.22 

9.5 

1600 

0.0684 

0.1467 

0.47 

2.23 

9.8 

400 

0.0037 

9.6 

warn 

125-50-30-2 























TLL 

30 

0.58 

0.0605 

0.4197 

0.14 

1.12 

14.1 

0.0815 

0.4330 

0.19 

1.22 

16.4 

1650 

0.0987 

0.4483 

0.22 

1.31 

18.8 

1050 

0.0172 

17.6 


BLT 

30 

0.58 

0.0600 

0.4178 

0.14 

1.12 

14.0 

0.0860 

0.4317 

0.20 

1.25 

16.7 

1650 

0.1134 

0.4665 

0.24 

1.38 

21.3 

1050 

0.0274 

19.0 

Bsfl 

TRL 

30 

0.58 

0.0605 

0.1357 

0.45 

2.14 

9.4 

0.0658 

0. 1425 

0.46 

2.21 

9.6 

1650 

0.0704 

0.1529 

0.46 

2.21 

10.0 

1050 

0.0046 

9.8 

9pff^9 

BRT 

30 

0.58 


0.1360 

0.49 

2.36 

9.3 

0.0660 

0.1425 

0.46 

2.22 

9.6 

1650 

0.0702 

0. 1546 

0.45 

2.18 

10.1 

1050 

0.0042 

9.8 


125-50-40-1 

TRL 

40 

0.77 

0.0590 

0.4155 

0.14 

1.06 

19.2 

0.0782 

0.4355 

0.18 

1.15 

22.4 

140 

0.1170* 

0.4960 

0.24 

1.30 

35.2 

630 

0.0388 

28.8 

61.6 

TLT 

LRL 

40 

40 

0.77 , 
0.77 : 

0.0560 

0.0590 

0.4046 

0.1210 

0.14 

0.49 

1.05 

2.28 

18.6 

11.9 

No Grow 
0.0637 

th 

0.1328 

0.48 

2.24 

12.5 

140 

140 

0.0792 

0.0660 

0.4130 

0.1384 

0.19 

0.48 

1.18 

2.23 

21.9 

12.8 

630 

630 

0.0232 

0.0023 

20.3 

12.7 

>36.9 

3.7 

LLT 

40 

0.77 

0.0606 

0.1255 

0.48 

2.26 

12.1 

0 . 0640 

0.1310 

0.49 

2.29 

12.4 

140 

0.0670 

0.1385 

0.48 

2.26 

12.8 

630 

0.0030 

12.6 

4.8 

125-50-40-2 























TLT 

31 

0.77 

rarc! 7 ’»i 

0.4125 

0.15 

1.07 

14.5 

0.0680 

0.4125 

0.16 

1.11 

20.2 

1345 

0.0798 

0.4205 

0.19 

1.17 

16.8 

1000 

0.0118 

16.2 

11.8 

BLL 

40 

0.77 

0.0612 

0.4180 

0.15 

1.07 

19.3 

0.0810 

0.4300 

0.19 

1.17 

22.1 

1345 

0.0990 

0.4504 

0.22 

1.25 

25.3 

1000 

0.0180 

23.7 

18.0 

TRT 

31 

0.77 

0.0564 

0.1360 

0.41 

1.95 

9.85 

0.0580 

0.1392 

0.41 

1.96 

12.8 

1345 

0.0620 

0.1475 

0.42 

1.97 

10.2 

1000 

0.0040 

10.1 ! 

4.0 

BRL 

40 

0.77 

0.0594 

0.1375 

0.43 

2.02 

12.7 

0.0620 

0.1430 

0.43 

2.03 

13.0 

1345 

0.0698 

! 0.1575 

0.44 

2.07 

13.8 

1000 

0.0078 

i3.4 ; 

7.8 

125-75-30-1 

TLT 

30 

0.58 

0.0910 


0.15 

1.14 

18.6 

0.0945 

0.0604 

0.16 

1.15 

19.2 

200 

0.1320* 

0.6205 

0.21 

1.29 

28.7 

1000 

0.0375 

23.9 

37.5 

BLL 

30 

0.58 

0.0930 

0.6041 

0.15 

1.14 

18.7 

0.0970 

0.6094 

0.16 

1.15 

19.3 

200 

0.1300* 

0.6314 

0.21 

1.27 

27.1 

1000 

0.0330 

23.2 

33.0 

TRT 

30 

0.58 

0.0921 


0.48 

2.32 

11.5 

0.0930 

0.1917 

0.49 

2.32 

11.6 

200 

0.0970 

0.1992 

0.49 

; 2.33 

12.0 

1000 

0.0040 

11.8 

4.0 

BRL 

30 

0.58 

0.0926 


0.45 

2.14 

12.2 

0.0945 

0.2098 

0.45 

2.16 

12.4 

200 

0.0985 

0.2195 

0.45 

2.15 

12.8 

1000 

0.0040 

12.6 

4.0 


tFirst two letters orient flaw in relation to sample number on panel (T - top, B = bottom, L = left, R ® right); last letter shows flaw 
orientation with respect to grain direction (T = transverse, L = longitudinal). 


-tFirst number represents material thickness in thousandths of an inch; secon d num ber represents nominal a/t ratio; third number indicates 
stress level in ksi. Minimum cyclic stress is zero, therefore, K « AK « FMa ^Tta/Q, where F = 1.0 + 0.12(l-a/2C) 2 , M is the back surface 
correction factor (Forman, R. G.: "Computer Analysis of Two-Dimensional Fatigue Flaw Growth Problem." NASA TMX-58086, Feb 1972); units 

of K are ksi^iiT. 

♦Indicates thru flaw. 


Table A-4 (oont) 
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vo 



























A-10 






EDM PREFLAW 




STAGE I GROWTH 




STAGE II GROWTH 



CYCLIC GROWTH 

VALUES 

Panelt 

Number 

a 

a/a 
/ ys 

a 

0 

2C 

0 

a/2C 

0 


V* 

0 

a l 

2c i 

a l /2C l 


V 

a 

a 2 

1 

a 2 /2C 2 

Q 2 

Ktv 

i 2 

N^or AN 

■ 

K +K "t 
i l i 2 

Aa 

AN 


(ksi) 


(in.) 

(in.) 




(in.) 

(in.) 





(in.) 

IPS 




(cycles) 


2 

y cycle / 

125-75-40-2 

TLT 

40 

0.77 ■ 

0.0955 

0.6260 

0.15 

1.08 

27.6 

0.0996 

0.6293 

0.16 

1.10 

28.6 

125 

0.1255* 

0.6438 

0.19 

1.19 

39.6 

125 

0.0259 

34.1 

207.2 

BLL 

40 

0.77 ' 

0 . 0900 

0.6010 

0.15 

1.08 

26.1 

0 . 1024 

0.6064 

0.17 

1.12 

29.8 

125 

0.1234* 

0.6300 

0.20 

1.19 

39.2 

75 

0.0210 

34.5 

280.0 

TRT 

40 

0.77 

0.0930 

0.1935 

0.48 

2.25 

16.2 

0.0937 

0 . 2000 

0.47 

2.19 

16.5 

125 

0.0960 

0.2043 

0.47 

2.19 

16.8 

125 

0.0023 

16.7 

18.4 

BRL 

40 

0.77 

0.0913 

0.1945 

0.47 

2.19 

16.2 

0.0934 

0.2015 

0.46 

2.17 

16.6 

125 

0.0980 

0.2060 

0.48 

2.22 

17.0 

125 

0 . 0046 

16.8 

36.8 


'First two letters orient flaw in relation to sample number on panel (T = top, B = bottom, L = left, R = right); last letter shows flaw 
orientation with respect to grain direction (T - transverse, L - longitudinal). 

'•'First number represents material thickness in thousandths of an inch; secon d num ber represents nominal a/t ratio; third number indicates 
stress level in ksi. Minimum cyclic stress is zero, therefore, K = AK. = FMo ^7Ta/Q, where F - 1.0 + 0 . 12(l-a/2C) 2 , M is the back surface 
correction fa ctor (Forman, R. G.: "Computer Analysis of Two-Dimensional Fatigue Flaw Growth Problem.” NASA TMX-58086, Feb 1972); units 

of K are ksi ^in. 

^Indicates thru flaw. 


Table A-4 (ooncl) 














TI-V 





EDM PREFLAW 



STAGE I 

GROWTH 



STAGE II GROWTH 

CYCLIC GROWTH VALUES 

Panelt 

Number 

0 

of a 
/ ys 

a 

0 

2C 

0 

a/2C 

0 

H 

K C 

a l 

2c i 

m 


K ft 
1 1 

yi 

a 2 

2C 2 

V 2C 2 

H 

K ft 
L 2 

N 2 orAN 


HEH 

Aa 

AN 











■ 



1 




■ 



■ 

■ ■ 

(P in ‘ ) 


(ksi) 


(in.) 

(in.) 




(in.) 

(in.) 





(in.) 

(in.) 




(cycles) 



\ cycle J 

PB-WL-25-FL 























15W 

15 

0.86 

0.0295 

0.1931 

0.15 

1.05 

4.9 

0.0295 

0.1931 

0.15 

1.05 

4.9 

7000 

0.0295 

0.1931 

0.15 

1.05 

4.9 

3000 

H,G. 

4.9 


50W 

15 

0.86 

0.0280 

0.0654 

0.43 

1.98 

3.3 

0.0280 

0.0654 

0.43 

1.98 

3.3 

7000 * 

0.0280 

0.0654 

0.43 

1.98 

3.3 

3000 

N.G, 

3.3 


15H 

15 

0.52 

0.0282 

0.1945 

0.14 

1.13 

4.6 

0:0282 

0.1945 

0.14 

1.13 

4.6 

7000 

0.0282 

0.1945 

0.14 

1.13 

4.6 

3000 

N,G, 

4.6 


50H 

15 

0.52 

0.0280 

0.0640 

0.44 

2.12 

3.2 ‘ 

0.0280 

0.0640 

0.44 

2.12 

3.2 

7000 

0.0280 

0.0640 

0.44 

2.12 

3.2 

3000 

N.G, 

3.2 


PB-WL-25-FT 










> 













15W 

30 

0.92 

0.0296 

0.1990 

0.15 

1.02 

9.9 

0.0325 

0.1995 

0.16 

1.05 

10.2 

1800 

0.0423 

0.2073 

0.20 

1.15 

11.1 

1300 

0.0098 

10.6 

7.5 

50W 

30. 

0.92 

0.0292 

0.0654 

0.45 

2.03 

6.6 

0.0307 

0.7000 

0.44 

2.00 

6.8 

1800 

0.0325 

0.0725 

0.44 

2.04 

7.0 

1300 

0.0018 

6.9 

1.4 

15H 

30 

0.58 

0.0286 

0.1960 

0.15 

1.12 

9.3 

0.0309 

0.1960 

0.16 

1.15 

9.5 

1800 

0.0357 

0.2012 

0.18 

1.20 

10.0 

1300 

0.0048 

9.8 

3.7 

50H 

30 

0.58 

0.0298 

0.0640 

0.47 

2.23 

6.4 

0.0310 

0.0660 

0.47 

2.25 

6.5 

1800 

0.0332 

0.0704 

0.47 

2.26 

6.7 

1300 

0.0022 

6.6 

1.7 

PB-WL-50-FL 























15W 

15 

0.86 

0.0563 

0.4025 

0.14 

1.02 

7.1 

0.0575 

0.4025 

0.14 

1.03 

7.2 

1250 

0.0594 

0.4025 

0.15 

1.04 

7.3 

2750 

0.0019 

7.2 

0.7 

50W 

15 

0.86 

0.0578 

0.1310 

0.44 

2.03 

4.7 

0.0578 

0.1310 

0.44 

2.03 

4.7 

1250 

0.0581 

0.1310 

0.44 

2.04 

4.7 

2750 

0.0003 

4.7 

0.1 

15H 

15 

0.52 

0.0605 

0.4195 

0.14 

1.13 

7.8 

0.0605 

0.4195 

0.15 

1.13 

7.8 

1250 

0.0605 

0.4195 

0.14 

1.13 

7.8 

2750 


7.8 


50H 

15 

0.52 

0.0575 

0.1276 

0.45 

2.17 

4.6 

0.0575 

0.1276 

0.45 

2.17 

4.6 

1250 

0.0582 

0.1276 

0.46 

2.20 

4.6 

2750 

0.0070 

4.6 

0.3 

PB-WL-50-FT 























15W 

30 

0.92 

0.0581 

0.4050 

0.14 

1.01 

14.6 

0.0593 

0.4200 

0.14 

1.00 

14.9 

300 

0.0672 

0.4200 

0.16 

1.04 

15.8 

210 

0.0079 

15.3 

37.6 

50W 

30 

0.92 

0.0595 

0.1320 

0.45 

2.05 

9.5 

0.0595 

0.1340 

0.44 

2.02 

9.6 

300 

0.0600 

0.1380 

0.43 

1.98 

9.8 

210 

0.0005 

9.7 

2.4 

15H 

30 

0.58 

0.0566 

0.4280 

0.13 

1.10 

14.8 

0.0730 

0.4520 

0.16 

1.16 

19.0 

300 

0.0865 

0.5340 

0.16 

1.16 

26.3 

210 

0.0135 

22.7 

64.3 

50H 

30 

0.58 

0.0587 

0. 1300 

0.45 

2.16 

9.2 

0.0595 

0.1350 

0.44 

2.12 

9.4 

300 

0.0612 

0.1495 

0.41 

1.98 

9.9 

210 

0.0017 

9.6 

8.1 


**Four flaws were placed in each panel, WL indicated weld was always oriented parallel to nominal 30 ksi loading direction, next two numbers, 
represent a/t ratio, FL means flaw oriented parallel to loading direction and FT means flaw is transverse to loading; flaw numbers show 
nominal aspect ratio of 0.15 or 0.50 and flaw location H for heat affected zone, W for weld zone. 


+Note yield strength values were: 

17.5 ksi for weld zone, transverse to load 

32.5 ksi for weld zone, parallel to load 

28.7 ksi for heat affected zone, transverse to load 

52.0 ksi for heat affected zone, parallel to load (same as parent metal) 


Table AS Flaw Growth Data for 0.125-in. Welded 2219-T87 
Aluminum - 2:1 Biaxial Loading 






















A-12 


Panel 

(flaw*) 

Number 

a 

(ksi) 


EDM PREFLAW ! 

STAGE I GROWTH 


STAGE II GR0 

WTH 


CYCLIC 

GROWTH V 

ALUES 

a/a 
/ ys 

a 

o 

(in.) 

2C 

o 

(in. ) 

a/2C 

Q o 

K i + 

o 

a i 

(in. ) 

2c i 
(in. ) 

a i/ 2C i 

<4 


N i 

(cycles) 

a 2 

(in. ) 

2C 2 

(in.) 

a 2/ 2C 2 

«2 


N 2 or AN 
(cycles) 

V a i 

(in.) 

K +K 7 
2 

Aa 

AN 

gN 

1S-060-50-30 

-4 

28 

0. S4 

0.0295 

0.0667 

0.44 

2.12 

6.0 

0.0297 

0.0667 

0.45 

2.14 

6.0 

3524 

0.0328 

0.0693 

0.47 

2.27 

6.1 

876 

0.0031 

6.1 

3.5 

5 

28 

0. 54 

0.0316 

0.0670 

0.47 

2.26 

6.0 

0.0320 

0.0680 

0.47 

2.25 

6.0 

3524 

0.0340 

0.0685 

0.50 

2.38 

6 . 1 

876 

0.0020 

6 . 1 

2 . 3 

7 

31 

0 60 

0.0304 

0.0681 

0.45 

2.14 

6.2 

0.0315 

0.0711 

0.44 

2.13 

6.8 

3524 

0.0322 

0.0728 

0.44 

2.12 

7.0 

876 

0.0007 

6.9 

0 . 8 

-8 

31 

0 60 

0.0323 

0.0695 

0.46 

2.23 

6.8 

0.0323 

0.0695 

0.46 

2.23 

6.8 

3524 

0.0323 

0.0695 

0.46 

2.23 

6 . 8 

876 

0.0000 

6 . 8 

0. 0 

-9 

36 

0 69 

0.0305 

0.0650 

0.47 

2.25 

8.1 

0.0323 

0.0685 

0.46 

2.21 

8.3 

3524 

0.0338 

0.0700 

0.56 

2.14 

8.3 

876 

0.0015 

8.3 

1.7 

-10 

36 

0.69 

0.0305 

0.0650 

0.47 

2.25 

8.1 

0.0325 

0.0693 

0.47 

2.25 

8.3 

3524 

0.0350 

0.0712 

0.49 

2.36 

8. 3 

876 

0.0025 

8.3 

2 . 9 

2S-060-50-30 

-4 

26 

0* 50 

0.0310 

0.1913 

0.16 

1.16 

9.0 

0.0332 

0.1925 

0.17 

1.18 

9.2 

2650 

0.0370 

0.1938 

0.19 

1.23 

9.5 

950 

0.0038 

9.4 

4.0 

5 

?6 

0.50 

0.0350 

0.1953 

0.18 

1.20 

8.4 

0.0356 

0.1953 

0.19 

1.23 

8.5 

2650 

0.0386 

0.1953 

0.20 

1.25 

8.7 

950 

0.0030 

8.6 

3. 2 

-7 

28 

0.54 

0.0296 

0.2190 

0.14 

1.10 

8.8 

0.0327 

0.2190 

0.15 

1.13 

9.2 

2650 

0.0406 

0.2190 

0.19 

1.22 

9.8 

950 

0.0079 

9.5 

8 . 3 

8 

?9 

0. 56 

0.0290 

0.2231 

0.13 

1.09 

9.4 

0.0318 

0.2255 

0.14 

1.11 

9.7 

2650 

0.0382 

0.2276 

0.17 

1.17 

10.3 

950 

0 . 0064 

10.0 

6 . 7 

-9 

31 

0.60 

0.0294 

0.2052 

0.14 

1.12 

9.3 

0.0330 

0.2090 

0.16 

1.15 

9.7 

2650 

0.0450 

0-2110 

0.21 

1.29 

10.6 

950 

0.0120 

10.1 

12.6 

-10 

35 

0.67 

0.0297 

0.2041 

0.15 

1.12 

12.9 

0.0343 

0.2085 

0.16 

1.16 

10.6 

2650 

0.0531 

0.2185 

0.24 

1.37 

14.7 

950 

0.0188 

14.0 

19 . 8 

3S-060-50-30 

-4 

23 

0.44 

0.0272 

0.0630 

0.43 

2.08 

4.7 

0.0272 

0.0630 

0.43 

2.08 

4. 7 

3000 

0.0285 

0.0640 

0.45 

2.14 

4.8 

2300 

0.0013 

4.8 

0.6 

5 

23 

0.44 

0.0305 

0.0642 

0.48 

2.28 

4.8 

0.0305 

0.0642 

0.48 

2.28 

4.8 

3000 

0.0323 

0.0650 

0.50 

2.38 

4. 8 

2300 

0.0018 

4 . 8 

0 . 8 

-7 

28 

0.54 

0.0561 

0.1405 

0.40 

1.94 

8.8 

0.0587 

0.1440 

0.41 

1.97 

9.0 

3000 

0.0668 

0.1513 

0.44 

2.12 

9.2 

2300 

0.0081 

9.1 

3.5 

-8 

28 

0.54 

0.0291 

0.0689 

0.42 

2.04 

6.2 

0.0312 

0.0710 

0.44 

2.11 

6.3 

3000 

0.0335 

0.0717 

0.47 

2.2 

6.2 

2300 

0.0023 

6.3 

1 . 0 

-9 

35 

0.67 

0.0290 

0.0700 

0.41 

2.00 

8.0 

0.0305 

0.0732 

0.42 

2.01 

8.1 

3000 

0.0344 

0.0805 

0.43 

2.1 

8. 4 

2300 

0.0039 

8. 3 

1. 7 

-10 

35 

0.67 

0.0526 

0.1427 

0.37 

1.81 

11.7 

0.0546 

0.1473 

0.37 

1.82 

11.9 

3000 

0.0721 

0.1773 

0.41 

2.0 

12.4 

2300 

0.0175 

12.3 

7 . 6 

4S-125-50-30 

-4 

20 

0.38 

0.0600 

0.1332 

0.45 

2.16 

8.6 

0.0600 

0.1367 

0.44 

2.11 

8.7 

2000 

0.0600 

0.1709 

0.35 

1.74 

8.8 

1500 

N.G. 

8.8 

N.G. 

_5 

20 

0.38 

0.0600 

0.1103 

0.54 

2.20 

7.6 

0.0600 

0.1103 

0.54 

2.20 

7.8 

2000 

0.0600 

0.1547 

0.39 

1.89 

8.0 

1500 

N . G , 

8.0 

N.G. 

-7 

29 

0.56 

0.0644 

0.1435 

0.45 

2.15 

9.0 

0.0650 

0.1445 

0.45 

2.16 

9.1 

2000 

0.0792 

0.1566 

0.51 

2.37 

9.4 

1500 

0.0142 

9.2 

9.5 

-8 

29 

0.56 

0.0590 

0.1395 

0.42 

2.04 

9.0 

0.0600 

0.1414 

0.42 

2.04 

9.1 

2000 

0.0660 

0.1571 

0.42 

2.03 

9.4 

1500 

0.0060 

9.2 

4.0 

-9 

32 

0.62 

0.0545 

0.1421 

0.38 

1.87 

10.4 

0.0565 

0.1482 

0.38 

1.86 

10.6 

2000 

0.0769 

0.1801 

0.43 

2.06 

11.6 

1500 

0.0204 

11.1 

13.6 

-10 

32 

0.62 

0.0582 

0.1388 

0; 42 

2.02 

10.3 

0.0610 

0.1428 

0.43 

2.06 

10.4 

2000 

0.0742 

0.1902 

0.39 

1.90 

11.9 

1500 

0.0132 

11.2 

8.8 

5S-125-50-30 

-4 

20 

0.38 

0.0630 

0.1300 

0.48 

2.32 

5.9 

0.0635 

0.1318 

0.48 

2.31 

6.0 

3700 

0.0645 

0.1330 

0.48 

2.32 

6.0 

828 

0.0010 

6.0 

1.2 

-5 

20 

0.38 

0.0600 

0.1306 

0.46 

2.20 

5.9 

0.0610 

0.1315 

0.46 

2.22 

6.0 

3700 

0.0635 

0.1352 

0.47 

2.25 

6.0 

828 

0.. 0025 

6.0 

3.0 

-7 

27 

0.52 

0.0580 

0.1456 

0.40 

1.93 

8.5 

0.0594 

0.1500 

0.40 

1.92 

8.7 

3700 

0.0645 

0.1516 

0.43 

2.05 

8.7 

828 

0.0051 

8.7 

6.2 

-8 

28 

0.54 

0.0574 

0.1430 

0.40 

1.94 

8.5 

0.0632 

0.1525 

0.41 

2.00 

8.8 

3700 

0.0754 

0.1-95 

0.42 

2.02 

9.7 

828 

0.0122 

9.3 

14. 7 

-9 

32 

0.62 

0.0530 

0.1450 

0.37 

1.80 

10.5 

0.0595 

0.1470 

0.40 

2.00 

10.6 

3700 

0.0684 

0.1500 

0.46 

2.20 

10.7 

828 

0.0089 

10.6 

10.7 

-10 

32 

0.62 

0.0560 

0.1434 

0.39 

1.90 

10.3 

0.0680 

0.1655 

0.41 

2.00 

11.1 

3700 

0.0790 

0.1800 

0.44 

2.10 

11.6 

828 

0.0110 

11.4 

13.3 

6S-125-50-30 

-4 

21 

0.40 

0.0650 

0.3760 

0.17 

1.18 

9.0 

0.0694 

0.3800 

0.18 

1.20 

9.3 

2200 

0.0791 

0.3917 

0.20 

10.0 

14.5 

2000 

0.0097 

9.7 

4.9 

-5 

21 

0.40 

0.0617 

0.3061 

0.20 

1.26 

8.3 

0.0624 

0.3061 

0.20 

1.27 

8.6 

2200 

0.0658 

0.3061 

0.21 

9.2 

13.0 2000 

0.0034 

8.8 

1.7 

-7 

26 

0.50 

0.0710 

0.4462 

0.16 

1.15 

12.0 

0.0735 

0.4462 

0.16 

1.17 

12.1 

2200 

0.0982 

0.4462 

0.22 

13.2 

15.9 

2000 

0.0247 

12.8 

12.4 

-8 

26 

0.50 

0.0602 

0.4456 

0.14 

1.09 

11.4 

0.0620 

0.4456 

0.14 

1.11 

11.5 

2200 

0.0890 

0.4456 

0.20 

13.5 

15.5 2000 

0.0270 

12.4 

13.5 

-9 

32 

0.62 

0.0580 

0.4332 

0.13 

1.09 

14.3 

0.0599 

0.4385 

0.14 

1.10 

14.5 

2200 

0.0964 

0.4490 

0.21 

16.7 

15.8 

2000 

0.0365 

15.6 

18.3 

-10 

26 

0.50 

0.0588 

0.4268 

0.14 

1.10 

11.5 

0.0605 

0.4302 

0.14 

1.11 

11.6 

2200 

0.1098 

0.4516 

0.24 

14.4 

16.4 

2000 

0.0493 

13.0 

24 . 6 


*Flaw number relates to flaw location in Figure V-ll and V-12 . 

iMinimum cyclic stress is zero; therefore K - MC - FHo where F = 1.0 + 0.12(l-e/2C) 2 , M is the back surface correction factor (Forman - 

NASA TM X-58086): units of K are ksl-^in. 


Table A- 6 Tabular Listing of Flaw Growth Bata for Biaxially 
Loaded Stiffened Panels ( Uniform Stress Field) 




Panel (Flaw*) 
Number 

a 

0 

(in.) 

2C 

0 

(in.) 

a l 

(in.) 

2c i 

(in.) 

a 2 

(in.) 

2C 2 

(in.) 

Aa/AN 

(p-in . /cycle) 

AKt 

(ksi ^in . ) 

1 

1S-060-50-30 










-1 

0.0298 

0.0680 

N.G.it 

N.G. 

N.G. 

N.G. 

— 

— 

— 

-2 

0.0317 

0.0674 

0.0426 

0.1019 

0.0464 

0,1125 

4.3 

9.0 

32.4 

-3 

0.0320 

0.0705 

N.G . 

N.G. 

N.G. 

N.G. 

— 

— 

— 

-6 

0.0280 

0.0695 

N.G . 

N.G. 

0.0285 

0.0704 

>0.6 

— 

— 

2S-060-50-30 










-1 

0.0330 

0.1450 

0.0358 

0.1460 

0.0460 

0.1480 

10.7 

13.6 

42.5 

-2 

0.0316 

0.1584 

0.0395 

0.1800 

0.0583 

0.2060 

19.7 

17.3 

47.2 

-3 

0.0338 

0.1995 

0.0354 

0.2000 

0.0382 

0.2050 

2.9 

7.4 

22.4 

-6 

0.0360 

0.2050 

0.0370 

0.2060 

0.0400 

0.2060 

3.1 

7.6 

22.6 

3S-060-50-30 










-1 

0.0320 

0.0635 

N.G. 

N.G. 

N.G. 

N.G. 

— 

— 

— 

-2 

0.0342 

0.0635 

0.0373 

0.0680 

0.0521 

0.0880 

6 . 4 

10.8 

40.6 

-3 

0.0324 

0.0658 

0.0325 

0.0658 

0.0350 

0.0680 

1.1 

5,5 

25.4 

-6 

0.0437 

0.0628 

N.G. 

N.G. 

N.G. 

N.G. 

— 

— 

— 

4S-125-50-30 










-1 

0.0660 

0.1350 

N.G. 

N.G. 

N.G. 

N.G. 

— 

— 

— 

-2 

0.0585 

0.1310 

0.0690 

0.1500 

0.0785 

0.2010 

6.3 

10.8 

29.3 

-3 

0.0616 

0.1355 

0.0624 

0.1365 

0.0665 

0.1370 

2,7 

7.1 

23.1 

-6 

0.0630 

0.1370 

N.G. 

N.G. 

0,0641 

0.1370 

>0.7 

— . 

— 

53-125-50-30 










-1 

0.0631 

0.1312 

0.0638 

0.1335 

0.0711 

0.1382 

8,8 

12,5 . 

40,3 

-2 

0.0611 

0.1330 

0.0680 

0.1590 

0,0815 

0.1900 

16.3 

16.2 

43.7 

-3 

0.0606 

0.1373 

0.0626 

0.1395 

0.0681 

0.1430 

6.6 

11.0 

34.7 

-6 

0.0580 

0.1358 

0.0588 

0.1364 

0.0603 

0.1375 

1,9 

6.0 

19.5 

63-125-50-30 

-1 

0.0692 

0.4165 

0.0778 . 

0.4165 

0.0963 

0.4165 

9.2 

12,8 

25.6 

-2 

0.0605 

0.4020 

0.0635 

! 0.4020 

0.1130 

0.4020 

24.7 

18,8 

37.0 

-3 

0.0786 

0.4148 

0,0800 

0.4148 

0.0970 

0.4148 

8.5 

12.3 

24.6 

-6 

0.0600 

0.4120 

0.0733 

0.4120 

‘ 0.0845 

0,4120’ 

5.6 

10.2 

21.2 

*Flaw numbers reflect location 

as shown 

in Figs 

V-ll and V-12, 





+AK from average biaxial growth rate curve and using obtained Aa/AN. 
**a computed from AK = FM o ^ira/Q . 
ttNo Growth, 


Table A- 7 Tabular Listing of Flaw Growth Data for Stiffened Panels } 
Panels } - Nonuniform Stress Field 
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1.0 INTRODUCTION 


The purpose of the photoelastic coating analysis of the aluminum 
panels was oriented to: 

1. Establish and define the stress distribution in the center 
bays of six aluminum panels. The stress distribution on 
both the front and rear sides was analyzed. 

2. Determine the level of stress concentration in the corners 
and fillet radii adjacent to the center bay. 

3. Determine the stress distribution in the stiffeners. 
Determination of the above information was accomplished on each 
of the six panels for three loading conditions: 

1 . 1:0 Uniaxial 1 oadi ng 

2 . 1:1 Biaxial loading 

3. 2:1 Biaxial loading I 

Application of plastic and photoelastic coating measurements 
were conducted by Photolastic, Incorporated. Martin-Marietta 
constructed the test facility and applied the loads per their 
specifications. 

2.0 INSTRUMENTATION 

2 . 1 Plastic Application 

Type PL-1* photoelastic coating was applied to the entire 
central bay of the six test panels defined on Martin dwgs. 
163-72-027 sheets 1 and 2. In general, particular attention 
was given to the fillet and corner areas of the panels. 
Figure 1 shows the area which was coated. Preparation of 
the contour sheets was in accordance with Photolastic, 
Incorporated's bulletin IB-P-310 (Instructions for Molding 
and Contouring Photo elastic Sheets, See Appendix). 

# /' . J 
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*Photol asti c , Inc., Malvern, Pa. 



Calibration of the individual plastic sheets established 
the photoelastic strain sensitivity. Bonding of the plasti 
to the panels was accomplished using type PC-1 reflective 

- - adhesive*. Figure 2 shows pertinent dimensions for each 

of the six test panels. 

2.2 Pol ar i scope 

Photoelastic measurements were made using an 031 Reflective 
Polariscope and the Oblique Incidence Attachment*. Photo- 
graphic recordings of the resulting strain patterns were 

made with a Nikon camera. 

1 . 

3.0 TEST PROCEDURE 

i 

The panels were installed in the loading fixture as shown in 
Figures 3A and 3B. Once a panel was installed, it was loaded in 
Uniaxial Tension (Vert.: Horiz.:: 1:0) 

Biaxial Tension (Vert.: Horiz.:: 1:1) 

Bijaxial Tension (Vert.: Horiz. : : 2:1) 
the chronological sequence of testing was: 

Spec No . Dates Tested 

5 11/14/72 - 11/15/72 

6 11/16/72 

5 4 11/17/72 

1 2 11/17/72 - 11/18/72 

v ‘ 1 1 1/18/72 - 1 1 /20/72 

3 11/21/72 

The first test panel (#5) was subjected to uniaxial loads of 
130, 86.6 and 43.3 kips. Experimental data showed that the 
specimen's response to external load was essentially linear. 

' ' I , 

*Photol asti c , Inc., Malvern, Pa. B ~ 3 



The following table summarizes the kips load when photoelastic 
measurements were made: 



LC 

AD CONDITION (KIP$) 


PANEL 

1 : 0 

1 : 1 

2:1 

1 

V = 65 & 25 
H = 0 & 0 

V = 65 & 25 
H = V 

V = 65 & 25 1 

H = 1/2 V 

2 

V = 65 & 25 
H = 0 & 0 

V = 65 & 25 
H = V 

V = 65 & 25 
H = 1/2 V 

3 

V = 65 & 25 
H = 0 & 0 

V = 65 & 25 
H = V 

V = 65 & 25 
H = 1/2 V 

4 

V = 130 & 43.3 
H = 0 & 0 

V = 1 00 & 43.3 
H = V 

V = 1 30 & 43.3 
H = 1/2 V 

5 

V = 130, 86.6 & 0 
H = 0, 0, & 0 

V = 1 30, 86.6 & 43.3 
H = V 

V = 130 & 43.3 
H = 1/2 V 

6 

V = 130 & 43.3 
H = 0 & 0 

V = 130 & 43.3 
H = V 

V = 130 & 65 
H = 1/2 V 


The load induced photoelastic signal (aN) was obtained as the 
difference between measurements (N) taken at the tabulated upper 
and lower load levels. For example, on panels 1 through 3 (1:0 
loading) : 

AN = N 65k " N 25k 

Here AN represents the signal induced by a 40 kip load. Martin's 
interest was to determine stress values induced by 130 kip on 
specimens 4, 5, and 6 and 91.6 kips on specimens 1, 2, and 3. 

All experimentally determined stresses were simply extrapolated 
11 nearly from the actual test load induced values to stress values 
corresponding with the expressed 130 kip and 91.6 kip levels. 

Data was collected at specific points designated by the coordinate 
system of Figure 4. The origin (HO, VO) was defined by the inter- 
section of the horizontal and vertical ribs. Points of measure- 
ment ( H2 , V3 for example) were described in terms of inches from 
the horizontal rib axis and vertical rib axis (that is 2 in. in 
a direction parallel to the horizontal rib and 3 in. parallel to 
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the vertical rib). Further definition of points of measurement 
are illustrated by the lower sketch of Figure 4. Lines identified 
as "A", "B", "C M and "D" define: 

"A" the top edge of the ribs 

"B" the tangent line between ribs and fillet 

"C" the center (45°) of the fillet 

"D" the tangent line between the central bay and the fillet 
For example, a location 1 inch from the origin (HO, VO) along the 
top edge of the vertical rib is designated as A1 . Further, point 
B 1/2 indicates a location on a point one half inch removed from 
the corner along the "B" line. Similar reasoning is used for 
lines "C" and "D". 




where "f" is established by calibration and is a function of 
coating thickness at the points of measurement. In thin skin 
applications, one finds that it becomes necessary to apply a 
correction factor (C-|) which will account for the reinforcement 
effect influenced by the plastic coating. This correction was 
taken into account and all tabulated data presented in this re- 
port are in final form. Determination of separate values of 
principal stresses (cr-| andc^) can be accomplished with an addi- 
tional photoelastic reading. This is accomplished by utilizing 
an oblique incidence attachment which is readily attachable to 
the reflecting polariscope. Thus, the emerging Tight from the 
polar i scope is passed obliquely through the plastic establishing 
an additional oblique incidence fringe order (N 0 ). In this 
investigation the oblique measurements were used only to check 
the degree of biaxiability of the stress field. 

5.0 TEST RESULTS 

Before entering a discussion of the resulting test data, one 
should become familiar with the coordinate system defined on 
Figure 4. Due to the symmetry of the problem a detailed analysis 
was conducted on one quadrant of the central bay. 

5 . 1 Uniaxial Loading (1:0 Loading) 

Figures 5A through 10B show the photoelastic coating pattern 
resulting from uniaxial loading. Tables IA, IB, and IC show 
the experimental results for the required vertical loads. 

5 . 2 Biaxial Loading (1:1 Loading) 

The photoelastic coating patterns are shown in Figures 1 1 A 
through 16B and the experimental results are shown in Tables 
1 1 A , 1 1 B , and IIC. 



5 . 3 Biaxial Loading (2:1 Loading) 

Figures 17A through 22B show the photoelastic coating 
patterns and Tables 1 1 1 A , 1 1 1 6 , and 1 1 1 C show the experi- 
mental results. 

6.0 DISCUSSION OF TEST RESULTS 

6.1 Uniaxial Loading (Tables IA, through IC & Figs. 5A through 10B) 
The test results, tabulated in Tables IA through IC, demon- 
strate negligible bending through the thickness of the 
central panel bays. The measured stress values at H2V2, 

H2V6 , H6V2 , and H6V6, on the front and back surfaces, show 
good agreement (approximately + 1 ksi or ± 3% of the nominal 
applied tensile stress). Further, the experimental data 
showed negligible stress gradients between the panel front 
face and the outside edge (top) of the vertical rib. This 
was particularly true in the central portion of the vertical 
rib (locations A4 through D4, A3 through D3, and A2 through 
D2 ) . Stress gradients near the corner, formed by the inter- 
secting vertical and horizontal ribs, were very steep. The 
experimental results show the stress levels at AO, BO, and 
CO to be measurably less than the nominal vertical rib 
stresses at mid-span. Reductions of 25% to 60% were observed. 
Stresses at DO were comparable to the nominal values of 
tensile stress in the panel faces and vertical ribs. The 
horizontal ribs experienced a ' 1 ow compressive stress along 
lines A and B. Measurements' were made along lines A and B 
on specimen 3 only. The experimentally determined stresses 
along line A (-10 ksi to -13 ksi) compared favorably with 
the predictable value -10 ksi. Stresses along line D 
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adjacent to the horizontal rib were comparable with the 
nominal tensile stresses in the central bay. The horizontal 
ribs generated a disturbance in the stress field on the back 
face (see Figures 5 through 10) with a maximum value occurring 
at location H0V0 (approximately 30% higher than the nominal 
tensile stress). Stresses along the VO line (H 1/2 VO 
through H4V0) were 10 to 25% greater than the nominal tensile 
stress values. The data indicate that the nominal vertical 
panel stress (specimen to specimen) ranged between a high 
of 35 ksi in specimen 1 to 29 ksi in specimens 3 and 6. The 
following observations are in order: 

-The corner stresses (corner formed by intersection 
ribs) were lower than the nominal tensile stress through 
the panel . 

-Bending stress through the panel thickness was 
generally less than + 3% of the applied membrane 
stress. 

-Back face stresses (immediately behind the horizontal 
ribs) were 10 to 25% greater than the nominal values 
measured in the central portion of the panel. 

-The maximum experimental stress was observed on the 
back face at coordinate H0V0 and was approximately 
30% greater than the nominal tensile stress. 

-Oblique incidence measurements established the stress 
field at H0V0, H1/2V0, H1V0, and H4V0 to be uniaxial. 

6.2 Equal Biaxial Loading (Tables 1 1 A to IIC, Figs. 1 1 A to 16B) 
Normal incidence measurements provide an immediate evalua- 
tion of equal biaxial tension in the central panel. Deviation 
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from the zero fringe (black) indicates lack of symmetry 
between vertical (° v ) and horizontal (<^) stress. Dif- 
ferences between o v and were generally less than 4 ksi 
or approximately 10% of the nominally applied stress. As 
in the uniaxial loadings, the stress gradients in the corner 
of the ribs were very steep. The data»in general, do not 
indicate a stress concentration in this area. Location 
CO of specimen 3 proved to be an exception to this general 
observation. No clear and obvious reason for this exception 
is postulated; however, it can be attributed to the general 
nonsymmetri cal stress distribution, present in varying degrees 
of all specimens. Stresses along line A of the ribs (mid- 
span) were generally in the 20 to 25 ksi range. This is in 

2 

qualitative agreement with the predictable value of 3 the 
nominal membrane stress in the panel. The photoelastic 
patterns on the back face provided a simple indication of 
the nonsymmetrical stress system existing in the specimens. 

A pure biaxial load would necessarily produce a symmetrical 
strain distribution. None of the specimens exhibited a 
truly symmetrical pattern; however, the pattern of Specimen 
3 (Figure 13B) approaches symmetry while the pattern on the 
back of Specimen 6 (Figure 16B) is a dramatic example of 
lack of symmetry. This 1:1 loading condition proved to be 
the most difficult to repeat. Further evidence of non- 
symmetrical loading was indicated by observed changes in the 
isoclinic angle (direction of principal stress) which 
occurred between the 25 kip and 65 kip load levels. In 
any event, the following observations are in order: 
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-The corner stresses (corner formed by intersecting 
ribs) were lower in magnitude than the nominal mem- 
brane stresses in the central section of the ribs. 

-The ribs produced stress risers on the back face. 
Stresses perpendicular to the rib directions were 
generally 1 5. to 30% greater than the nominal panel 
membrane stress. 

6.3 Two to One Biaxial Load (a v : ah:: 2:1) 

(Tables 1 1 1 A through IIIB, Figures 17A through 22B) 
Measurements on the front and back faces (points H2V2, 
H2V6, H6V2, and H6V6) again confirmed bending through 
the panel thickness to be low (approximately - 1.5 ksi). 
Stresses in the central region of the vertical fillet 
and rib were essentially uniform with an indication of 
reduced stress levels at locations AO, BO, and CO. The 
experimental stresses at DO were generally comparable to 
the membrane stresses in the panel faces and the strain 
pattern indicated no stress concentration at the location. 
Assuming that a v = 2 ah, it is possible to predict the 
vertical and horizontal rib stresses as .84 o v and .17 a v 
respectively*. The rib stresses along lines A & B are 
seen to be in general agreement with these predictions. 
There was a unique behavior along the fillet of the hori- 
zontal rib which was not observed during the uniaxial and 
equal biaxial loadings. The principal stress directions 
were observed to be measurably different from vertical 
and horizontal (tending towards 45°) along line C» 
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*E = 10(10)® and v = 


.33 for aluminum 




This behavior is reasonable and is illustrated here, 
a-j in the horizontal rib is directed along the length of 
the rib while a-j (o v as applied) is vertical in the central 
bay of the panel. Clearly then, the a-j stress trajectory 
must experience a 90° rotation as it crosses through the 
fillet region and is illustrated above by the dotted line. 
This unique situation was not present in the other loading 
cases. As in the previous loading conditions the back 
face stresses were highest behind the horizontal rib along 
line VO with the vertical stresses 15 to 30% greater than 
the nominal applied membrane stress. 

The following observations are in order: 

-Bending stress through the central portion of the 
panel was approximately - 1 1/2 ksi. 

-The corner stresses (corner formed by intersecting 
ribs) at AO, BO and CO were lower in magnitude 
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than the nominal vertical membrane stress. 

-The horizontal ribs produced stress risers on the back 
face (Line VO). Stresses perpendicular to the hori- 
zontal ribs were 15 to 30% greater than the nominal 
vertical membrane stress. 

7.0 GENERAL OBSERVATIONS AND SUGGESTIONS 

A. Analysis of the experimental findings indicates that the 
ribs create high stresses on the back face of the panel. 

The stress increase (up to 30%) occurs directly behind the 
ribs in a direction perpendicular to the rib length. 

B. The results show that the stress field in the central section 
of the panel can be treated as plane stress. The bending 
component was generally less than - 2 ksi for all loadings. 

C. The stress gradients in the corners formed by intersecting 
ribs were very steep. In any case, the general conclusion 
derived from consideration of all tests is that the stresses 
at AO, BO, and CO are less than the nominally applied mem- 
brane stress. 

D. The above observations are based on measurements made in one 
gradient of the specimen. Lack of loading symmetry most 
certainly produced different stress fields in each quadrant. 
However, observation of the overall photoelastic patterns 
suggests that the above are generally applicable to the 
enti re panel . 

E. An overall review of all test results leads to the observation 
that the fillet radii, rather than rib thickness, was the 
dominant factor influencing localized stress disturbances. 
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F. It is suggested that the above observations be confirmed 
by further strain gage studies. Figure 23 indicates loca- 
tions where strain gage measurements would provide useful 
comparative data. 
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2 - Panel Dimensions 




FIGURE 3A PANEL IN TEST FIXTURE 



FIGURE 3B OVERALL VIEW OF TEST FIXTURE 
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FIGURE 5B -Panel 1, UNIXIAL LOAD 
REAR SURFACE 





FIGURE 6A-PANEL 2 UNIAXIAL LOAD 
FRONT SURFACE 



FIGURE 6B-PANEL 2 UNIAXIAL LOAD 
REAR SURFACE 
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FIGURE 7B PANEL 3 UNIAXIAL LOAD 
REAR SURFACE 


B-29 





FIGURE 8A- PANEL 4 UNIAXIAL LOAD 
FRONT SURFACE 



FIGURE 8B-PANEL 4 UNIAXIAL LOAD 
REAR SURFACE 
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FIGURE 9A-PANEL 5 UNIAXIAL LOAD 
FRONT SURFACE 



FIGURE 9B -PANEL 5 UNIAXIAL LOAD 
REAR SURFACE 
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FIGURE 10-A PANEL 6 UNIAXIAL LOAD 
FRONT SURFACE 



FIGURE 10-B PANEL 6 UNIAXIAL LOAD 
•REAR SURFACE 
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FIGURE 1 1 B- PANE L 1 -BIAXIAL LOAD (1:1) 
REAR SURFACE 
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FIGURE 12A- PANEL 2 BIAXIAL LOAD (1:1) 
FRONT SURFACE 



FIGURE 12B- PANEL 2 BIAXIAL LOAD (1:1) 
REAR SURFACE 


B-34 







FIGURE 1 3A- PANEL 3 BIAXIAL LOAD (1:1) 
FRONT SURFACE 
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FIGURE 14A- PANEL 4 BIAXIAL LOAD (1:1) 
FRONT SURFACE 



FIGURE 1 4 B -PANE L 4 BIAXIAL LOAD (1:1) 
REAR SURFACE 
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FIGURE 15A 
PANEL 5 

BIAXIAL LOAD (1:1) 
FRONT SURFACE 




FIGURE 15B 
PANEL 5 ' 

BIAXIAL LOAD (1:1) 
CORNER VIEW 



FIGURE 15C 
PANEL 5 

BIAXIAL LOAD (1:1) 
REAR VIEW 
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FIGURE IGA- PANEL 6 BIAXIAL LOAD (1:1) 
FRONT SURFACE 



FIGURE 1 6B- PANEL 6 BIAXIAL LOAD (1:1) 
REAR SURFACE 
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FIGURE 17A PANEL 1 BIAXIAL LOAD (2:1) 
FRONT SURFACE 



FIGURE 17B PANEL 1 BIAXIAL LOAD (2:1) 
REAR SURFACE 
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FIGURE 18B PANEL 2 BIAXIAL LOAD (2:1) 
REAR SURFACE 


FIGURE 18A- PANEL 2 BIAXIAL LOAD (2:1) 
FRONT SURFACE 
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FIGURE 1 98- PANEL 3 BIAXIAL LOAD (2:1) 
REAR SURFACE 


B-41 


FIGURE 20A- PANEL 4 BIAXIAL LOAD (2:1) 
FRONT SURFACE 



FIGURE 20B- PANEL 4 BIAXIAL LOAD (2:1) 
REAR SURFACE 



FIGURE 21A -PANEL 5 BIAXIAL LOAD (2:1) 
FRONT SURFACE 
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FIGURE 22B-PANEL 6 BIAXIAL LOAD (2-1) 
REAR SURFACE 
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